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ABSTRACT 


Since  i960  there  has  been  a  great  deal  of  interest 
in  experiments  involving  electron  tunneling  into  super¬ 
conductors.  This  new  technique,  which  is  very  simple  in 
principle,  provides  a  direct  method  of  observing  the 
energy  gap  in  the  density  of  electron  states  of  a  super¬ 
conductor.  More  recently  this  technique  has  also  been 
applied  to  the  study  of  detailed  variation  in  the  density 
of  states  in  a  superconductor.  The  theory  of  electron 
tunneling,  as  it  applies  to  the  study  of  superconductors, 
is  discussed. 

In  order  to  be  able  to  observe  small  variations 
in  the  density  of  states  of  superconductors  a  new  method 
for  measuring  the  dynamic  conductance  of  tunnel  junctions 
had  to  be  devised.  This  method  involves  a  very  sensitive 
bridge  technique  and  is  described  in  detail.  Since  some 
of  the  experiments  of  interest  had  to  be  performed  at 
temperatures  below,  1°K  a  He^  cryostat  was  used.  A 
description  of  the  design  and  operation  as  well  as  proposed 
modifications  of  this  apparatus  is  given.  Although  the 
cryostat  has  so  far  only  been  used  in  tunneling  experiments, 
it  was  designed  in  such  a  way  that  it  would  be  used  to 
investigate  other  phenomena,  particularly  transport  proper¬ 
ties  in  metals. 
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Initial  experiments  involved  measurements  of 
the  energy  gap  in  superconducting  lead.  During  these 
experiments  we  gained  familiarity  with  tunneling  technique 
and  refined  our  experimental  method.  The  energy  gaps  of 
superconducting  zinc,  aluminum,  and  indium  were  also 
measured.  All  of  these  are  found  to  be  in  good  agree¬ 
ment  with  the  theory  of  Bardeen,  Cooper,  and  Schrieffer 
(1957b),  and  with  the  exception  of  zinc  for  which  no 
published  tunneling  data  were  found,  with  the  values 
obtained  by  other  workers  using  tunneling  techniques. 

The  discovery  of  structure  in  the  density  of 
electron  states  of  superconducting  lead  by  C-iaever,  Hart, 
and  Megerle  (1962)  and  by  Rowell,  Chynoweth,  and  Phillips 
(1962),  stimulated  us  to  search  for  similar  effects  in 
other  metals.  Investigation  by  Adler  and  Rogers  (1963) 
of  the  density  of  states  of  superconducting  indium  showed 
that  if  such  structure  is  present  it  must  be  smaller  by 
more  than  an  order  of  magnitude  than  that  observed  in 
lead.  However  new  and  thus  far  unexplained  structure  in 
the  density  of  states  of  indium  was  found. 

Finally,  a  number  of  tunneling  experiments  in¬ 
volving  superconductors,  and  which  should  lead  to  more 
definitive  results,  are  proposed.  A  new  method  for  inves¬ 
tigation  of  the  Fermi  surface  of  normal  metals  using 
electron  tunneling  is  also  suggested. 
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CHAPTER  I 

THEORETICAL  APPROACH  TO  SUPERCONDUCTIVITY 


A.  Historical  Introduction 

Superconductivity  was  discovered  by  Onnes  (1911) 
while  investigating  the  temperature  variation  of  the 
electrical  resistance  of  mercury  wires  at  liquid  helium 
temperatures.  He  found  that  the  resistance  of  the  wires 
he  was  using  became  residual  (i.e.  reached  a  constant 
value  determined  by  the  scattering  of  electrons  by  lattice 
imperfections  or  impurities)  at  about  l4°K  and  then  quite 
abruptly  at  about  4°K  the  resistance  dropped  to  zero. 

Since  that  time  the  complete  disappearance  of  resistance 
has  been  observed  in  about  two  dozen  pure  metals  and  over 
one  hundred  alloys.  Shortly  after  the  discovery  of  this 
superconductive  state  it  was  found  that  by  applying  a 
magnetic  field  greater  than  a  certain  critical  value,  H  , 
resistive  behaviour  could  be  restored.  It  was  further 
found  that  the  critical  field  increased  as  the  temperature 
was  lowered,  from  an  initial  value  of  zero  at  the  super¬ 
conducting  transition  temperature,  T  ,  to  a  finite  value 
at  the  absolute  zero.  This  behaviour,  along  with  the 
almost  discontinuous  change  in  the  specific  heat  at  the 
transition  temperature  indicated  that  the  change  which 
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a  metal  undergoes  In  becoming  a  superconductor  is  a  phase 
transition  in  the  thermodynamic  sense.  This  transform¬ 
ation  is  of  the  second  order  if  there  is  no  magnetic  field 
present,,  while  a  first  order  transition  takes  place  in 
non-zero  magnetic  fields. 

Application  of  Maxwell's  electromagnetic  equations 
to  materials  of  infinite  conductivity  indicates  that  the 
magnetic  flux  density.,  B,  is  a  constant  in  such  media.  If 
this  constant  were  other  than  zero  the  normal  to  supercon¬ 
ducting  transition  would  be  dependent  on  the  previous 
history  of  the  specimen.,  and  reversible  thermodynamics 
could  not  then  be  applied  to  this  problem.  However  the 
experiments  of  Meissner  and  Oschenfeld  (1933)  showed  that 
indeed  the  superconducting  state  is  one  of  perfect  dia¬ 
magnetism  (i.e.  B  =  const.  =  0)  and  hence  the  transition 
is  a  thermodynamically  reversible  one.  In  the  period 
between  1933  and  1956  a  number  of  "two  fluid  models" 
flourished.  These  models  treated  the  electrons  in  the 
superconductor  as  if  they  were  constituents  of  two  inter¬ 
penetrating  gases.  One  of  these  gases,  the  normal  com¬ 
ponent,  consisted  of  ordinary  electrons  not  unlike  those 
considered  in  the  Sommerfeld  model  of  a  metal.  The 
second  component  consisted  of  "super-electrons",  which  did 
not  contribute  to  the  entropy  of  the  electron  gas  and  were 
responsible  for  the  superconducting  properties.  The 
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theories  of  the  Londons  ( cf .  London  (1950))  gave  a  fairly 
good  macroscopic  explanation  of  the  electromagnetic  pro¬ 
perties  of  superconductors  "based  on  a  modification  and 
expansion  of  Maxwell's  equations.  Many  phenomenological 
two  fluid  models  describing  the  thermodynamic  behaviour 
of  superconductors  appeared  starting  with  the  Casimir  and 
Gorter  (193^-)  model.  All  of  these  have  since  been  shown 
to  be  special  cases  of  a  more  general  model  due  to  Marcus 
(1952  a,  1952b,  1953)  and  Lewis  (1956). 

Very  little  progress  toward  a  microscopic  theory 
was  made  before  1950.  Until  that  time  two  problems  in 
particular  eluded  explanation:  (i)  only  the  electronic 
properties  and  not  the  lattice  properties  were  affected  by 
the  transition,  (ii)  many  metals  in  spite  of  their  vastly 
different  band  structure  became  superconductors. 

B.  Attempts  at  a  Microscopic  Theory  of  Superconductivity 

Significant  progress  towards  a  microscopic 
theory  of  superconductivity  occurred  with  the  almost 
simultaneous  discovery  by  Reynolds  et  al .  (1950)  and 
prediction  by  Frihlich  (1950  a,b)  of  the  isotope  effect. 

This  effect  demonstrated  that  in  many  superconductors 
the  superconducting  transition  temperature,  T  ,  is  depend¬ 
ent  on  the  isotopic  mass,  M,  according  to  the  relation 

T  M^/^ 
c 


const . 


(1.1) 


. 
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for  a  given  element.  Thus  it  was  conclusively  shown  that 
although  the  lattice  remained  unaffected  by  the  transition 
it  does  play  an  important  role  in  superconductivity. 

Frbbfich  (1950  a,b)  proposed  a  theory  of  super¬ 
conductivity  ba&ed  on  an  attractive  electron-electron 
interaction.  The  mechanism  of  this  interaction  was  due  to 
the  exchange  of  phonons'  between  the  interacting  electrons. 
This  interaction  is  possible  because  the  ions  are  neither 
rigidly  bound  in  the  lattice,  nor  infinitely  heavy.  Thus 
an  electron  passing  through  the  lattice  slightly  distorts 
the  periodic  potential.  Another  electron  traveling 
through  the  distorted  (polarized)  lattice  does  not  exper¬ 
ience  the  same  periodic  potential  which  the  first  electron 
saw;  in  this  way  the  two  electrons  interact  via  the  lattice. 
FrtSlich  further  showed  that  if  the  energy  difference 
between  the  electron  states  in  such  phonon  exchange  inter¬ 
actions  is  less  than  the  phonon  energy  the  interaction  is 
attractive.  He  was  also  able  to  establish  the  criterion 
for  the  occurrence  of  superconductivity,  i.e.  that  if 
such  an  attractive  electron-electron  interaction  dominates 
over  repulsive  Coulomb  interaction  between  electrons,  a 
metal  becomes  superconducting.  He  was  able  to  explain. 


* 


Strictly  speaking  "virtual"  phonons  since  superconduc¬ 
tivity  presumably  exists  at  T  =  0°K. 
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qualitatively,  why  poor  conductors  (having  strong  electron- 
phonon  coupling)  became  superconductors  while  good  con¬ 
ductors  (e.g.  group  one  metals)  did  not.  His  theory  also 
accounted  foe  the  isotope  effect.  The  quantitative  features 
of  the  theory  failed  however  because  these  involved  calcu¬ 
lating  the  energy  of  the  superconducting  state  using  pertur¬ 
bation  theory.  The  difficulty  of  this  approach  is  evident 
if  we  consider  that  the  energies  of  the  normal  and  super¬ 
conducting  states  are  both  about  one  electron-volt  per 

atom,  while  the  energy  difference  between  them  is  less 

-7 

than  10  of  this  value.  The  calculation  of  the  energy 
difference  by  perturbation  methods  could  not  be  carried 
out  with  sufficient  accuracy.  Attempts  at  a  similar 
theory  by  Bardeen  (1950  a,b)  using  the  variational  approach 
also  resulted  in  mathematical  difficulties.  Various 
attempts  by  other  authors  to  produce  a  microscopic  theory 
of  superconductivity  are  reviewed  by  Bardeen  (1956). 

C,  Experimental  Evidence  of  an  Energy  Gap  in  Superconductors 
In  the  ensuing  seven  years  a  great  deal  of  experi¬ 
mental  evidence  accumulated  indicating  an  energy  gap,  2A, 
in  the  density  of  electron  states  of  a  superconductor 
such  that  no  excited  electronic  states  exist  whose  energy 
is  less  than  2A  above  the  ground  state. 


Si  5 
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(i)  Corak  et  al.  (1956)  observed  that  the  low 
temperature  specific  heat  of  superconductors  is  proportional 
to  exp ( -A/kT) . 

(ii)  Blevins,  Gordy,  and  Fairbank ( 1955 )  observed 
that  the  absorption  of  electromagnetic  radiation  increases 
sharply  when  the  frequency.,  co,  is  greater  than  a  critical 
value  2A/fi.  They  also  demonstrated  that  the  energy  gap, 

2A,  increases  from  zero  at  the  transition  temperature  to 

a  limiting  value  of  2A  £  3-5  kTc,  at  absolute  zero. 

(iii)  Experiments  by  Glover  and  Tinkham  (1956, 
1957)  involving  infrared  transmission  in  thin  superconduct¬ 
ing  films  also  indicated  an  energy  gap. 

(iv)  Mackinnon  (1955)*  Mason  and  Bommel  (1955)* 
and  Morse,  Tamerkin  and  Bohm  (1955)  all  showed  that  ultra¬ 
sonic  attenuation  decreases  sharply  in  the  superconducting 
state . 

Many  other  experiments  gave  further  evidence 
for  an  energy  gap.  A  complete  and  lucid  account  of  these 
is  given  by  Biondi,  Forrester,  Garfunkel,  and  Satterthwaite 
(1958) . 

I).  The  BOS*  Theory 

This  theory,  although  in  some  respects  a  simp¬ 
lified  model,  is  able  to  account  for  most  of  the  phenomena 

*  Throughout  this  thesis  BCS  refers  to  the  theory  of 
Bardeen,  Cooper  and  Schrieffer  (1957  a,b) . 
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of  the  superconducting  state.  No  attempt  to  describe  it 
in  any  detail  will  be  made  here.,  instead  only  those  aspects 
pertinent  to  electron  tunneling  phenomena  will  be  dealt 
with.  We  begin  by  considering  the  behaviour  of  electrons 
in  a  normal  metal,  and  then  describe  how  this  differs  in 
a  superconductor  according  to  the  BCS  theory. 

It  is  indeed  striking  that  in  spite  of  the  complex 
band  structure  and  diversity  of  metals  which  become  super¬ 
conductors  there  is  great  similarity  in  their  behaviour 
below  the  transition  temperature.  We  know  from  Bloch's 
theorem  that  normal  metals  may  be  described  in  terms  of  one 
electron  wave  functions  of  the  form 


-ik  -r 


(1.2) 


where  k,s  are  the  four  quantum  numbers  corresponding  to 
wave  vector  and  spin  respectively,  and  Ugf  g  is  a  function 
having  the  periodicity  of  the  lattice.  Since  the  electrons 
are  Fermions,  the  many-electron  s  y  s  t  e  m  will  consist 
of  all  the  energy  levels  being  filled  up  to  the  Fermi 
level  (energy  and  wave  vector  kJV)  .  The  only  correlations 
between  electrons  dealt  with  in  the  Bloch  description  of 
a  metal  are  statistical  ones  arising  due  to  the  antisymmetry 
requirement  of  the  wave  function.  In  a  normal  metal  there 


We  deal  here  with  the  special  case  of  a  spherical 
Fermi  surface  since  that  is  the  only  case  treated  by 
BCS. 


' 
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will  be  a  large  number  of  single  particle  excitations 
(of  energy  %  kT)  in  which  an  electron  in  a  state 
will  be  replaced  by  one  in  a  state  k.  >  k^ .  Such 
excitations  may  be  conveniently  described  in  terms  of 
"quasi-particle"  energies,  i.e.  the  excitation  energy 
measured  from  the  Fermi  surface.  The  energy  difference 
between  the  ground  state  and  an  excited  state  is  then 
given  in  terms  of  the  sum  of  such  quasi-particle  energies, 
e(k).  Phenomena  such  as  the  electrical  resistance  of  a 
normal  metal  arise  through  interactions,,  in  particular 
the  electron-phonon  and  electron-impurity  interactions. 

Both  of  these  processes  provide  a  mechanism  by  which  the 
electrons  can  readily  be  scattered  from  one  momentum 
state  into  another.  Such  a  mechanism,  for  example.,  assures 
that  any  ordered  state,  such  as  a  current  carrying  state, 
will  be  quickly  randomized  in  the  absence  of  the  field 
which  produced  it.  Herein  lies  the  essential  difference 
between  the  normal  and  superconducting  state,  for  in  the 
superconducting  case  such  single  particle  excitations  are 
far  less  frequent  and  consequently  ordered  states  may 
persist  even  in  the  absence  of  the  fields  which  produced 
them.  This  difference  between  the  normal  and  supercon¬ 
ducting  states  is  assumed  to  be  due  to  electron  correlations. 
It  is  by  introducing  correlations  that  the  BCS  theory  is 
able  to  treat  superconductivity. 


,  ii)  . 
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Cooper  (1956)  showed  that  if  there  exists  a  net 
attractive  interaction  between  electrons,,  then  pairs  of 
electrons  may  form  a  bound  state  with  a  net  gain  in  energy 
over  the  Fermi  ground  state.  Such  binding  arises  from 
long  range  correlations  between  pairs  of  electrons.,  and 
is  strongest  between  pairs  of  opposite  spins  and  momenta 
(-k  |  ,  k 4  ;  where  the  arrows  indicate  spin  orientation). 
Unfortunately  as  pointed  out  by  Frblich  (1962)  no  simple 
physical  picture  of  these  correlations  appears  to  exist. 
Some  insight  into  the  physical  basis  for  such  binding  can 
however  be  obtained  through  an  argument  due  to  Frdlich 
(1962)  if  one  assumes  that  the  attractive  interaction 
exists  only  for  electrons  above  the  Fermi  surface.  Con¬ 
sider  two  electrons  that  in  the  absence  of  mutual  inter¬ 
action  would  be  above  the  Fermi  surface,  thus  having 
energies  slightly  above  .  Owing  to  the  attractive 
interaction  the  resultant  energy  of  the  pair  will  be 
somewhat  less  than  2e^.  The  pair  cannot  however  break 
up  into  two  single  particles  because  all  the  energetically 
accessible  states  are  occupied  (since  such  states  lie 
below  the  Fermi  surface) .  This  pairing  of  electrons 
accounts  for  the  different  features  of  the  normal  and 
superconducting  phases  since  it  takes  a  small  but  finite 
amount  of  energy  to  break  up  a  pair.  This  results  in  a 
high  degree  of  coherence  of  the 
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wave  functions  of  the  superconducting  state,  which  makes 
single  particle  exitations  difficult  in  the  superconducting 
state . 

The  BCS  theory  assumes  that  the  major  difference 
between  the  normal  and  superconducting  states  arises  due 
to  this  added  dynamical  correlation  between  such  electron 
pairs  with  opposite  spin  and  momentum.  They  make  no 
assumption  about  the  nature  of  the  interaction  giving 
rise  to  such  pairs,  save  that  there  is  some  constant 
attractive  interaction  at  temperatures  below  Tc  which 
vanishes  for  T  =  Tc .  On  this  basis  they  calculate  the 
free  energy  difference  between  the  normal  and  superconduct¬ 
ing  states.  Their  theory  results  in  an  energy  gap  2A 
centered  about  the  Fermi  surface.  This  energy  gap  is  given 
in  terms  of  a  non-linear  integral  equation  which  involves 
the  attractive  electron-electron  interaction,  -V,  the 
density  of  states  at  the  Fermi  surface  of  the  normal 
metal,  p(0),  and  an  average  phonon  frequency,  co  , 
which  they  choose  to  be  half  the  Debye  frequency.  They 
are  able  to  solve  this  equation  in  the  limit 

P ( o)  v  «  l  *  .  (1.3) 

The  energy  gap  2A  at  the  absolute  zero  Is  then  given  by 

2A  =  4  (fio))  e"1/P(°)V  (1.4) 

o  '  ' ave  ^  ' 


* 


This  is  the  weak  coupling  limit,  i.e.  where  the  electron- 
phonon  interaction  is  weak. 


m 


.. 
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and  the  single  particle  energy  with  respect  to  the  Fermi 
surface  is 

E  =  (e2  +  A2)1/2  .  (1.5) 

It  is  evident  that  (1.5)  reduces  to  the  appropriate  quasi¬ 
particle  energy.,  e,  for  a  normal  metal  since  in  that  case 
the  energy  gap  is  zero.  They  are  able  from  the  above  to 
deduce  the  thermodynamic  properties  of  superconductors, 
as  well  as  the  temperature  variation  of  the  energy  gap. 

They  obtain  the  relation  2Aq  =  3.52  kTc  which  is  in  reason¬ 
able  agreement  with  the  experimental  values.  Further 
extensions  of  this  theory  by  Mattis  and  Bardeen  (1958), 
and  by  Rickayzen  (1958,  1959)  are  able  to  account  for  the 
electromagnetic  properties  of  superconductors.  A  thorough 
review  of  the  theory  and  recent  experimental  results  is 
given  by  Bardeen  and  Schrieffer  (1961). 

It  should  be  stressed  that  since  the  nature  of 
the  attractive  interaction  does  not  explicitly  enter  in 
the  BCS  theory,  the  framework  of  the  theory  should  in 
principle  also  be  able  to  deal  with  the  problems  of  transi¬ 
tion  metal  superconductors  ( cf .  G-eballe,  Mathias,  Hull, 
and  Gorenzwit  1961)  in  which  the  attractive  potential  does 
not  arise  from  the  electron-phonon  interaction.  The  effect 
of  explicit  types  of  interactions  have  been  treated  by 
Swi h  art  ( 1962) . 
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E.  The  Density  of  Electron  States  In  a  Superconductor 


A  quantity  of  prime  interest  in  electron  tunnel 
ing  into  superconductors  is  the  density  of  states.  This 
may  he  obtained  directly  from  (1.5). 


dN(E) 

dE 

dN ( e )  d  e 
de  dE 

=  Pn(0)  ^ 

n  dE 

where  p  ( 0 )  is  the  density  of  states  of  the  normal  metal 
evaluated  at  the  Fermi  surface.  Carrying  out  the  differ¬ 
entiation  of  (1.5)  we  obtain 


PS(E) 

-  Pn(°)RP  -  W.6) 

V  E 

where  ps(e)  is  the  density  of  states  of  the  superconductor 
It  is  evident  from  (1.6)  that  the  density  of  states  of  the 
superconductor  is  singular  at  the  edges  of  the  energy  gap, 
where  E  =  A.  The  total  number  of  states  remains  unaltered 
by  the  transition  to  the  superconducting  state. 


RP  denotes  the  real  part 


' 
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CHAPTER  II 


A  SURVEY  OF  THE  THEORY  OF  ELECTRON  TUNNELING 
INTO  SUPERCONDUCTORS 


A.  Introduction 

Until  August  i960  far  infrared  transmission  and 
ultrasonic  attenuation  presented  the  most  direct  methods 
of  observing  the  superconducting  energy  gap.  At  that  time 
Giaever  (1960a)  published  the  results  of  his  first  experi¬ 
ments  on  electron  tunneling  into  superconductors.  His 
technique.,  simple  in  principle ^  provides  a  direct  method 
of  observing  the  energy  gap  of  superconductors.  Since  then 
this  technique  has  also  been  applied  to  determining  the 
density  of  electronic  states  in  superconductors. 

Basically  what  is  done  in  tunneling  experiments 
is  to  inject  electrons  of  various  energies  into  a  super¬ 
conductor;  due  to  the  energy  gap  those  electrons  having 
their  energies  in  the  gap  region  will  be  rejected.  Such 

'  O 

injection  is  done  through  a  very  thin  (<100  A)  insulating 
barrier  which  allows  one  to  vary  the  potential  difference 
between  the  metal  from  which  the  electrons  emanate  and  the 
one  into  which  they  are  injected.  In  this  way  the  injec¬ 
tion  energy  may  be  varied  at  will. 


' 
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Bt  A  Model  of  Tunnel ing* 

It  has  long  been  known  ( cf .  Sammerfeld  and 
Bethe  1933)  that  a  current  can  flow  through  a  thin  Insu¬ 
lating  film  between  two  metals  due  to  quantum-mechanical 
tunneling.  A  short  qualitative  description  of  the 
tunneling  process  is  given  below  for  the  following  three 
cases:  (i)  two  normal  metals,,  (ii)  a  normal  metal  and  a 
superconductor ,  and  (iii)  two  superconductors. 

(i)  Pig.  2.1  (a)  shows  two  normal  metals  sepa¬ 
rated  by  a  thin  insulating  film  at  T  =  0°K.  Energy  is 
represented  on  the  abscissa,,  while  the  density  of  states 
near  the  Fermi  level  is  shown  on  the  ordinate**,  the 
corresponding  current-voltage  (I-V)  characteristic  is 
also  shown.  The  transmission  probability  for  a  particle 
tunneling  through  a  potential  barrier  depends  exponentially 
on  both  the  square  root  of  the  barrier  height  and  the 
thickness  of  the  barrier.  If  the  applied  voltage  is  kept 
small,  neither  the  barrier  height  nor  its  thickness  are 
affected.  The  current  will  increase  linearly  with  applied 


*  The  treatment  given  .below,  in  the  main  parallels  that 
given  by  Glaever  (1961J  and  by  Shapiro,  Smith,  Nicol, 
Miles,  and  Strong ( 1962 ) .  Since  no  review  article  on 
the  subject  is  yet  available,  I  attempt  in  this  section 
to  summarize  the  theory  briefly,  stressing  some  of  the 
underlying  assumptions.  Some  of  the  integrations,  the 
results  of  which  are  merely  stated  in  the  current  lit¬ 
erature  are  carried  out  in  detail  in  Appendix  I. 

**  Only  the  energy  region  within  a  few  kT  of  the  Fermi 
surface  is  shown  in  Fig.  2.1  the  lower  levels  are 
assumed  to  be  completely  filled. 
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FIGURE  2.1 

Schematic  description  of  various  tunneling  processes. 
This  figure  shows  the  density  of  states.,  D(E),  near 
the  Fermi  surface  plotted  against  energy,  E.  The  cor¬ 
responding  current-voltage  curves  are  also  shown. 
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voltage  ,  since  the  number  of  electrons  which  flow  is 
proportional  to  the  potential  difference  between  the  two 

-V.  AZ. 

sides  of  the  barrier.  Since  kT  is  small  compared  with 
the  barrier  height  in  our  work  the  tunneling  process  is 
temperature  independent . 

(ii)  Consider  now  the  situation  which  would 
exist  at  T  =  0°K  with  one  of  the  metals  superconducting 
(Fig.  2.1(c)).  The  density  of  states  of  the  supercon¬ 
ductor  (metal  l)  has  an  energy  gap  2A^  centered  about  the 
Fermi  level.  Thus  no  current  will  flow  until  the  applied 
voltage  is  at  least  V  =  A-^/e.  Assuming  that  the  number  of 
electrons  tunneling  is  proportional  to  the  density  of  states., 
then  as  V  increases  (V>A^/e)  the  current  increases  rapidly 
at  first.  Finally.,  as  the  electrons  face  a  less  favourable 
density  of  states,  the  current  approaches  assymptotically 
that  obtained  when  both  metals  are  normal.  At  non-zero 
temperatures  (Fig.  2.1(d))  there  are  always  some  thermally 
excited  electrons  above  the  energy  gap.  These  will  give 
a  temperature  dependent  contribution  to  the  current  even 
for  applied  voltages,  V  <  A^/e . 


*  This  is  true  for  the  small  applied  oltages  involved. 
For  detailed  analysis  see  Holm  (1950* ** 

**  N.B.  For  zero  applied  voltage  the  Fermi  levels  in  both 
metals  lie  at  the  same  level.  For  an  applied  voltage, 
V,  the  Fermi  levels  are  displaced  by  an  amount  eV 
from  one  another,  e  being  the  electronic  charge. 
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(ill)  This  situation  is  one  in  which  both  metals 
are  superconducting  as  shown  in  Pig.  2.1(e)  and  (f ) .  Again 
we  consider  first  the  case  at  T  =  0°K.  As  the  applied 
voltage  is  increased  no  current  flows  until  V  =  (A^  +  A2)/e 
For  larger  applied  voltages  the  current  rises  sharply  at 
first.,  and  finally  approaches  the  I-V  curve  obtained  when 
both  metals  were  normal.  For  T  >  0°K  there  are  again  a 
few  normal  electrons  above  the  gap.  At  any  given  tempera¬ 
ture  their  relative  numbers  will  be  proportional  to  Ap/A2. 
These  normal  electrons  give  rise  to  a  current  flow  even 
for  the  smallest  applied  voltages,  this  current  increases 
until  V  =  (A-j  -  A2)/e.  As  the  voltage  is  further  increased 
only  the  same  number  of  electrons  can  flow,  but  these  now 
face  a  less  favourable  density  of  states,  and  the  current 
therefore  decreases  with  increasing  voltage.  This  part  of 
the  I-V  curve  is  called  the  "negative  resistance"  region. 
Finally  when  the  applied  voltage  reaches  a  value  of 
(Aj  +  A2)/e  the  current  again  rises  sharply  at  first  and 
then  gradually  approaches  the  I-V  curve  obtained  for  both 
metals  In  the  normal  state. 

We  shall  now  derive  the  I-V  characteristics  for 
tunnel  Junctions  based  on  the  model  outlined  above.  Tunnel 
ing  through  a  potential  barrier  is  treated  as  an  ordinary 
quantum  transition.  The  transition  probability  per  unit 
time  for  the  transition  from  an  occupied  state  on  one 
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side  of  the  barrier  to  an  unoccupied  state  on  the  other 
is  given  by: 

M2  f1(1-ff)pf  -  (2.1)* 

where  is  the  density  of  final  states.,  f^  is  the  probabil¬ 
ity  that  the  final  state  is  occupied.,  and  f^  is  the  prob¬ 
ability  of  occupation  of  the  initial  state.  M  is  the  matrix 
element  for  the  transition.  Before  proceeding  it  is  neces¬ 
sary  to  make  the  following  assumptions: 

(i)  M  is  independent  of  energy  for  the  small 
energy  differences  involved, 

(ii)  M  remains  unaltered  in  the  transition  from 
the  normal  to  the  superconducting  state . 

Bardeen  (196!)  has  examined  these  assumptions  from  the  many- 
body  point  of  view,  and  using  the  W.K.B.  approximation  has 
shown  them  to  be  reasonable.  A  more  detailed  theoretical 
justification  of  these  assumptions  has  been  given  by  Cohen, 
Palicov,  and  Phillips  (1962).  It  should  be  noted  however 
that  Harrison  (1961)  has  also  shown  that  in  normal  metals 
these  assumptions  can  not  be  justified  on  the  basis  of 
the  independent  particle  model  alone.  If  these  two  assump¬ 
tions  are  made,  then  the  only  variable  that  enters  in  the 
transition  probability  is  the  density  of  final  states. 


*  See  any  standard  quantum  mechanics  text,  eg.:  L.I.  Schiff. 
Quantum  Mechanics,  McGraw-Hill  Co.,  New  York  (1955); 

Sec.  29. 
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If  we  now  consider  a  voltage,  V,  applied  to  one 
of  the  metals,  equation  (2.1)  enables  us  to  calculate  the 
tunneling  current.  The  net  tunneling  current  will  be 
given  by  the  sum  of  the  two  one  way  currents,  i.e. 

1=  (47re)/h/°°  |  M  |2  (p  (E+eV)  f(E)  p  (E)  [  1-f  (E+eV)  ]?  dE 
-00  ^  ' 

°°  Of 

~(4t re)/fij  |  M  f  |  p2(E+eV)  f(E+eV)  p^E)  [  1-f  (E)  ]jdE.  (2.2) 

The  subscripts  1  and  2  refer  to  metals  1  and  2  respectively, 

E  is  the  energy  measured  with  respect  to  the  Fermi  level, 
and  e  is  the  electronic  charge.  The  extra  factor  of  two 
is  present  due  to  the  two  possible  spin  states  of  the  elec¬ 
trons.  The  assumptions  about  M  allow  us  to  take  it  out  of 
the  integral  and  then 

00 

I  =  a  /  p„(E+eV)  p  (E)  [f (E)-f (E+eV) ]dE,  (2.3) 

-00 

where  a  =  (47re|  M  p)/h . 

For  the  case  when  both  metals  are  normal  and  at  T  =  0°K 
this  reduces  to: 

Jnn  =  aPl(0)  P2(°)eV^  (2.4)* 

where  p1(0)  and  p2(0)  are  the  densities  of  states  of 
metals  1  and  2  at  the  Fermi  surface.  This  expression  is 
valid  only  for  small  applied  voltages  as  is  the  case  in 


* 


See  Appendix  1(a)  for  Integration  (section  (i)) 


, 
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our  work.  A  further  assumption,  that  the  density  of  states 
in  the  vicinity  of  the  Fermi  surface  is  constant,  had  to 
be  made  in  order  to  obtain  equation  2.4. 

Similarly  for  the  case  in  which  one  metal  is  in 
the  normal  state  while  the  other  is  in  the  superconducting 
state  at  T  =  0°K  one  obtains:* ** 

Ins  =  0  for  eV  <  A  (2.5) 

Ins  =  ap1(0)p2(0)[  (eVU-A2]1/2  for  eV  >  A  (2.5a) 

G-iaever  and  Megerle  (1961)  have  evaluated  Ing  for  T  >  0°K 
and  eV  < A  based  on  the  BCS  theory.  They  obtain: 

Ins  =  ap1(0)p2(0)  2A  Z(-l)m+1  K1  ( mA/kT )  sinh(meV/kT) , 

m  (2.5b) 

where  is  the  first  order  modified  Bessel  function  of  the 
second  kind.  It  is  evident  from  the  above  that  the  super¬ 
conducting  energy  gap  as  well  as  its  variation  with  temper¬ 
ature  can  be  determined  by  tunneling  experiments.  The 

** 

results  of  such  experiments  are  in  good  agreement  with  the 
values  predicted  by  the  BUS  theory.  Good  qualitative 
agreement  is  also  obtained  between  the  experimental  I-V 


*  See  Appendix  1(a)  section  (ii)  for  derivation. 

**  cf.  Giaever  (1960a  and  1960b),  Giaever  and  Megerle  (1961), 
Giaever,  Hart,  and  Megerle  (1962),  Nlcol,  Shapiro,  and 
Smith  (i960),  Shapiro,  et  al.  (1962),  Townsend  and  Sutton 
(1962),  Taylor,  and  Bursteln  (1963). 


curves  and  those  predicted  using  the  BCS  density  of 
states  in  conjunction  with  the  model  of  tunneling  out¬ 
lined  above,  thus  giving  strong  support  to  the  BCS  theory. 
Having  made  the  assumption  that  the  matrix  element,  M, 
is  constant  the  current  becomes  proportional  to  the  density 
of  states.  Then  ratio  of  the  conductance,  (dl/dV)ns, 
when  one  of  the  metals  is  superconducting  to  the  conductance, 
(dl/dV)nn,  when  both  metals  are  normal  gives  the  relative 
change  in  the  density  of  states.  This  quantity,  called 
the  relative  conductance,  Gng  =  (dl/dV)  /(dl/dV)  ,  may 
be  calculated  from  equations  (2.4,  2.5.,  and  2.5a)^for 


T  =  OK: 


(dl/dV)  P  ?  wo 

Gns  =  7 - 7 — ~  =  (ev)[(ev)2  -  A2]'172.  (2.6) 

ns  dl/dV 

/  /  nn 


Experiments  of  Giaever  and  Megerle  (1961)  show  that  in 
spite  of  "kT  smearing"*  measured  values  of  relative  con¬ 
ductance  also  give  good  agreement  with  the  BCS  density  of 
states.  Thus  relative  conductance  measurements  on  tunnel 
junctions  consisting  of  a  normal  metal  and  a  superconductor 
provide  a  sensitive  tool  for  studying  the  energy  gap,  its 
temperature  dependence,  and  the  ratio  of  the  density 


The  physical  reason  for  this  smearing  being  that  for 
T  >  0°K  there  are  some  thermally  excited  normal  electrons 
above  the  gap. 
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of  states  of  the  superconducting  to  normal  phases  of  a 
superconductor.  In  the  temperature  region  below  1/2  Tc 
the  energy  gap  is  almost  temperature  independent  and  kT 
smearing  is  small  for  energies  eV  >  A.  Thus  equation 
(2.6)  may  be  used  for  comparison  with  experimental  data; 
this  avoids  laborious  numerical  integrations  on  a  computer 
which  are  required  to  calculate  the  relative  conductance 
for  T  >  0°K. 


We  now  turn  our  attention  to  the  case  of  two 


superconductors  separated  by  a  thin  insulating  film.  A 
scheme  by  which  numerical  evaluation  of  the  integrals 
arising  in  a  calculation  of  the  current  for  cases  in 
which  T  >  0°K  has  been  given  by  Shapiro,,  et  al.  (1962). 
For  simplicity  we  shall  limit  ourselves  here  only  to  the 
case  of  two  identical  superconductors  at  T  =  0°K.  The 
I-V  characteristic  for  this  case  is  given  below^for  a 
detailed  derivation  see  Appendix  I  (a)  section  (ill). 


I 


ss 


0  for  eV  >  2A 


(2.7) 


(2.7a) 


where  K  and  E  are  complete  elliptic  Integrals  of  the  first 
and  second  kind  respectively.  This  expression  has  been 


‘ 
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evaluated  on  the  IBM  1620  computer  and  a  table  of  relative 

conductance  as  a  function  of  voltage  is  given  in  Appendix 

II.  A  comparison  of  the  relative  conductance.  0  ,  for 

ss 

this  case  with  experiment  is  shown  in  Fig.  6.6. 

It  is  indeed  gratifying  that  good  qualitative 
agreement  is  obtained  between  this  simple  model  and  experi¬ 
ment.  Recalling  that  the  basic  assumption  of  the  model  was 
that  the  transition  probability  is  energy  independent , and  this 
led  to  the  current  being  proportional  to  the  densities  of 
states  of  the  superconductors.,  it  is  surprising  that  no 
density  of  states  variations  have  yet  been  detected  in 
normal  metals  using  tunneling  techniques.  Such  variations 
would  be  expected  to  be  rather  pronounced  in  some  of  the 
transition  metals.  Experimental  attempts  to  observe  such 
variations  (Giaever,  Hart,  and  Megerle  (1962))  in  nickel 
have  led  to  negative  results. 

C.  Modifications  to  the  BCS  Density  of  States  in  the 

Vicinity  of  the  Energy  Gap  as  indicated  by  Tunneling 
Measurements . 

With  one  exception  it  has  not  been  possible  to 
show  any  major  discrepancy  between  the  BCS  theory  and 
tunneling  experiments.  Equations  (2.7)  and  (2.7a)  show 
that  for  two  superconductors  at  T  =  0°K  there  ought  to  be 
a  current  jump  at  eV  =  2A.  Experiments  by  Giaever,  Hart, 
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and  Megerle  (1962)  on  various  tunnel  junctions  (particularly 
Sn-SnOx-Sn)  show  that  the  I-V  curve  has  a  finite  slope  at 
eV  =  2A.  This  slope  is  less  than  the  expected  after  the 
fact  that  the  measurements  are  made  at  T  >  0°K  is  taken  into 
account.  This  unexpected  slope  introduces  an  uncertainty 
in  the  determination  of  the  energy  gap  from  the  I-V  character¬ 
istics.  Giaever  et  al.  were  able  to  improve  the  agreement 
between  theory  and  experiment  by  assuming  that  the  energy 
levels  are  not  sharply  defined  but  instead  have  a  non- zero 
width.  Such  a  modification  tends  to  smear  the  density  of 
states  in  the  region  of  the  gap.  This  alteration  to  the 
BGS  theory  was  first  introduced  by  Hebei  and  Slichter  (1959) 
to  explain  the  results  of  nuclear  spin  relaxation  experi¬ 
ments  in  superconductors.  We  shall  now  describe  this 
modification  refering  to  it  as  HS  smearing  and  the  result¬ 
ing  density  of  states  as  the  HS  density  of  states  respec¬ 
tively.  The  change  is  introduced  into  the  theory  by 
replacing  the  sharp  energy  levels  of  BOS  by  others  having 
a  width  6.  6  is  later  adjusted  to  obtain  best  fit  with 

the  experimental  results.  The  energy  level  breadth  function 
used  has  a  width  26  and  a  height  1/26.  The  HS  density  of 
states  is  then  given  by: 

EH- 6 

PHS  =  25  ffecs  dE' 


(2.9) 


: 
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Carrying  out  the  integration  one  obtains: 


PHS  = 

0  for  | 

E  |  <  A  -  6 

(2.10) 

PHS  = 

1/2 6 [ (E+6)2 

-  A2]1/2  for  A  -  6  <  | 

E  |  <  A  +  6  (2.10a) 

PHS  = 

1/26  [[  (&h6)2 

-  A2 ] 1/2  _  [ ( E_  5 ) 2  _ 

a2]1/2} 

for 

|  E  |  >  A  +  6. 

(2.10b) 

It  can  be  seen  that  (2.10b)  goes  over  to  the  BCS  density  of 
states  as  the  energy,,  E<,  increases.  The  important  effect 
of  this  treatment  is  the  removal  of  the  singularity  in  the 
BCS  density  of  states.  Ctiaever  et  al.  found  improved 
agreement  between  theory  and  experiment  using  a  value  of 

6/A  =  .03.  If  the  smearing  is  attributed  to  finite  life- 

*  -11 
time  effects  ^  they  then  obtain  a  lower  limit  of  10  sec. 

for  the  lifetime. 

In  conclusion  it  is  evident  that  the  superconduct¬ 
ing  energy  gap  as  well  as  the  density  of  states  in  a  super¬ 
conductor  may  be  studied  by  the  technique  of  electron 
tunneling.  The  good  qualitative  agreement  between  such 
experiments  and  the  BCS  theory  is  one  of  the  major  triumphs 
of  that  theory.  The  discrepancy  which  does  arise  in  the 
vicinity  of  the  energy  gap  can  be  reduced  by  introduction 
of  a  level  breadth.  Further  discussion  is  given  in  the 
chapter  dealing  with  our  own  results. 


Assuming  the  lifetime  to  be  given  by  the  uncertainty 
principle ,  r  V6- 


D.  Multiphonon  Effects  Observed  in  Tunneling  Experiments 


C-iaever,  Hart,,  Megerle  (1962)  and  more  recently 
Rowells  Chynoweth,  and  Phillips  (1962)  have  reported  humps 
in  the  relative  conductance  curves  obtained  from  measure¬ 
ments  on  tunnel  junctions  in  which  at  least  one  of  the 
metals  is  lead.  Rowell.,  et  al .  were  able  to  resolve  har¬ 
monic  structure  in  the  conductance  curves  which  consisted 
of  peaks  occuring  at  energies: 

En  =  A  +  nEt  ,  (2.11) 

where  E^  is  the  energy  of  the  transverse  acoustic  phonon 
in  lead.,  n  is  an  integer,,  and  the  energy  En  =  eV  is  the 
electron  energy.  They  also  observe  that  the  magnitude 
of  the  peaks  decreases  with  increasing  voltage.  They 
attribute  this  observed  structure  to  multiphonon  effects 
(emission  or  absorption  of  one  or  more  phonons  by  the 
tunneling  electrons) .  In  arriving  at  their  explanation 
they  had  to  assume  an  Einstein  phonon  spectrum,  the  energy 
of  each  harmonic  oscillator  being  the  transverse  phonon 
energy.  Assuming  such  a  spectrum  a  simple  qualitative 
model  of  this  effect  can  be  given.  Consider  electrons 
tunneling  from  one  metal  to  another  and  having  some 
energy  E.  Even  if  there  are  many  available  states  of 
energy,  E,  only  those  electrons  having  the  appropriate 
spin  and  momentum  will  be  able  to  participate  in  the 


. 
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tunneling  process  and  thus  contribute  to  the  tunneling 
current.  Electrons  whose  energy  is  En  can  cross  the  barrier, 
emit  a  phonon,  and  thereby  reduce  their  energy  to  E1  =  A. 

In  this  manner  such  electrons  arrive  into  a  region  where 
the  density  of  states  available  to  them  is  high,  hence  by 
virtue  of  an  Einstein  spectrum  the  conduction  can  be  enhanced 
at  voltages,  V  =  E  /e.  The  reason  for  the  decreasing 
amplitude  of  successive  peaks  is  a  consequence  of  the 
fact  that  in  general  an  n  +  1  phonon  process  is  less  prob¬ 
able  than  an  n  phonon  process.  It  is  not  surprising 
that  such  structure  has  thus  far  only  been  observed  in 
lead  since  this  metal  has  the  strongest  electron-phonon 
interaction  of  all  superconductors. 


28. 


CHAPTER  III  • 

ELECTRICAL  APPARATUS  AND  MEASUREMENT  CIRCUITS 

A.  Introduction 

One  of  the  initially  appealing  aspects  about 
electron  tunneling  measurements  is  that  it  is  relatively 
simple  to  determine  the  I-V  curves  for  a  tunnel  junction. 
Since  I  is  of  the  order  of  microamperes ,  while  V  is  in  the 
millivolt  range;  a  good  quality  VTVM  will  suffice  to  measure 
both  I  and  V  (I  being  obtained  by  measuring  the  voltage 

developed  across  a  calibrated  resistor) .  Such  a  procedure 

Rf  R  r  1  •  1 

is  tedious.,  but  can  be  remedied  by  replacing  the  VTVM  with 
an  oscilloscope  and  simultaneously  displaying  V  on  the 
abscissa  and  voltage  proportional  to  I  on  the  ordinate. 

Two  ways  of  carrying  out  such  a  measurement  are  possible. 

The  first  of  these  involves  a  variable  dc  voltage  supply 
by  which  the  bias  across  the  tunneling  junction  may  be 
manually  varied.  The  second  method  uses  an  oscillator  to 
sweep  the  voltage  across  the  junction  and  displays  the 
I-V  curve  on  the  scope  screen  as  a  continuous  trace. 

Either  method  leads  to  accurate  results  provided  that 
four  terminal  measurements  are  made.  In  our  experiments 
we  have  employed  both  of  these  methods.  In  order  to 
record  this  data  we  employed  an  oscilloscope  camera  and 


. 
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negative  film.  The  use  of  negative  material  is  particu¬ 
larly  convenient  since  it  can  readily  be  enlarged  for 
closer  inspection.  When  using  the  oscillator  technique 
care  must  be  taken  to  keep  the  sweep  frequencies  low 
(below  a  few  cycles  per  second)  because  at  higher  frequencies 
junction  capacity  effects  become  visible. 

The  measurement  of  relative  conductance  on  the 
other  hand  is  not  quite  so  simple.  G-iaever,  Hart.,  and 
Megerle  (1962)  have  used  fairly  complex  circuitry  to 
record  the  dynamic  resistance  of  tunnel  junctions  directly 
on  an  X-Y  plotter.  They  carried  out  such  measurements  by 
superimposing  a  small  alternating  current  signal  on  the 
bias  current  through  the  junction  and  observing  the  ac 
voltage  developed  across  the  junction. 

We  have  used  a  completely  different  method  for 
obtaining  the  relative  conductance .  By  keeping  the 
apparatus  simple  we  have  been  able  to  obtain  much  higher 
sensitivity  than  other  workers  in  this  field,  although 
the  reduction  of  the  data  to  its  final  form  is  still 
somewhat  laborious. 

B,  The  Tunneling  Curve  Tracer 

A  complete  block  diagram  of  the  circuitry  is 
given  in  figure  3.1.  The  Hewlett  Packard  oscillator 
provided  a  triangular  waveform  which  was  used  for  the 


. 


r 


30. 


FIGURE  3.1 

Block  diagram  of  the  electrical  apparatus  used  in 
measuring  the  relative  conductance  and  the  current- 
voltage  curves  of  tunnel  junctions. 


(HEAVY  LINE  DENOTES  GROUND  WIRES) 
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automatic  sweep.  The  specimen  in  the  cryostat  was  connected 
to  the  external  apparatus  by  a  four  conductor  shielded 
cable.  Leads  1  and  4  are  the  current  leads;  2  and  3  are 
the  potential  leads.  The  heavier  black  lines  denote 
ground  bus  wiring.  Extreme  care  was  taken  to  ground  all 
equipment  and  shields  at  a  single  point  so  that  ac  pickup 
was  kept  to  a  minimum.  The  potentiometer  served  a  dual 
purpose.  With  the  switch  across  the  galvanometer  closed 
it  was  used  as  a  standard  voltage  source  to  calibrate 
the  oscilloscope  and  the  galvanometer  amplifier.  It  also 
served  as  a  bias  source  so  that  any  segment  of  the  I-V 
or  relative  conductance  curves  could  be  shifted  in  voltage 
for  a  more  detailed  display  on  the  oscilloscope  screen. 

For  relative  conductance  measurements  a  quantity *  <5V,  which 
is  the  voltage  difference  between  some  known  linear  I-V 
relation  and  the  specimen  voltage  at  the  same  current  was 
recorded  on  the  vertical  input.  Simultaneous  recording  of 
the  specimen  voltage,  V,  was  carried  out  on  the  horizontal 
input . 

The  galvanometer  amplifier,*  used  to  pre-amplify 
6V,  is  a  highly  stable  feedback  amplifier.  These  amplifiers 
are  discussed  extensively  by  Rogers  (1962)  but  the  one 


It  is  a  pleasure  to  acknowledge  my  gratitude  to  Mr. 

J.  S.  Rogers  for  designing  and  building  the  amplifier 
specifically  for  this  application. 
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used  here  is  somewhat  unconventional  in  that  the  output 
galvanometer  is  replaced  by  an  oscilloscope.  Excellent 
stability  is  achieved  due  to  high  feedback  and  the 

consequent  sacrifice  in  gain  is  not  serious  because  an 
amplification  factor  greater  than  10^  has  never  been 
required.  At  this  gain  a  voltage.,  6V,  of  one  microvolt 
produces  a  displacement  of  one  centimeter  on  the  oscillor- 
scope  screen. 

A  detailed  diagram  of  the  curve  tracer  appears 
in  Fig.  3.2.  As  can  be  seen  there  is  provision  for 
both  an  automatic  and  a  manual  sweep,  depending  on  the 
position  of  switch  Sw^ .  The  magnitude  of  the  voltage 
bias  across  the  specimen  can  be  controlled  by  adjusting 
SW2  and  the  250  ohm  variable  resistor.  The  ground  con¬ 
nections  are  all  made  at  terminal  8  and  the  1000  ohm 
potentiometer  can  be  adjusted  so  that  the  circuit  is  at 
some  level  above  ground  potential.  This  potentiometer  was 
adjusted  so  that  the  ac  pick-up  visible  on  the  oscillo¬ 
scope  screen  was  a  minimum.  The  I-V  curves  are  obtained 
by  the  four  terminal  method;  the  current  wiring  is  shown 
by  the  heavy  black  lines.  The  voltage  proportional  to  cur¬ 
rent  is  generated  across  R  ,  which  may  be  adjusted  to 
any  one  of  nine  values,  thus  providing  variable  current 
sensitivity. 


FIGURE  3.2 


Circuit  diagram  of  the  curve  tracer. 
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The  remainder  of  the  circuit  is  for  obtaining 
the  curves  of  6V  as  a  function  of  V,  from  which  the  relative 
conductance  may  be  calculated.  An  equivalent  circuit  for 
this  part  of  the  curve  tracer  appears  in  Fig.  3-3.  It  is 
essentially  a  bridge.,  one  arm  of  which  (the  tunnel  junc¬ 
tion)  is  a  non-linear  element.  If  we  let  Rg  be  the  input 
resistance  of  the  galvanometer  amplifier  then  we  may  write 
6V  =  i  R„ .  The  loop  equations  describing  the  circuit 

o  o 

may  then  be  written: 


e  -  V 


(3.1) 


e 


(3.2) 


o 


-(Ra+Rt>)±  +  Rb1o  +  (Ra+Rb+Rg);Lg 


(3.3) 


This  particular  set  of  loops  is  used  in  the  analysis 


since  the  form  of  V  as  a  function  of  i  is  not  known. 

Replacing  i  R  by  6V  and  combining  (3.1)  with  (3-2)  we 

§>  § 


have 


and 


(3.3a) 


Solving  these  equations  for  the  current  in  the  tunnel. 


junction  we  obtain 


' 


FIGURE  3.3 


Equivalent  circuit  of  the  relative  conductance 
measuring  part  of  the  curve  tracer. 
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i 


_1  /  Rb 
1  ^d 


(V+  6V)  +  [ 


(Ra+Rb) 


Rf 


+  1 


(3.5) 


Differentiating  (3-5)  with  respect  to  the  voltage  across 
the  junction  gives  the  dynamic  conductance, 


di 

dV 


R 


b 


■^d^a 


Rd(Ra+Rb)  Rd  d(5v) 
1  +  [  1  +  — — -  +  ^ 


Rb^h 


b 


dV 


(3.6) 


In  our  case  Ra  =  R^  =  R  so  (3.6)  becomes: 


di 

dV 


R 


d 


(i  +  [  +  j  a©)  . 

I  R  R  dV  ) 


(3.7) 


Note  that  Rd  in  Pig.  3.2  is  adjustable.  Let  us  designate 
the  setting  of  Rd  by  Rn  when  the  measurements  are  being 
carried  out  with  specimens  in  the  normal  state,  and  by 
Rs  for  the  measurements  with  superconducting  specimens. 

Let  <5v$  and  6Vn  denote  values  of  6V  obtained  from  measure¬ 
ments  with  superconducting  and  normal  specimens  respectively. 
Further  let  us  make  the  following  definitions: 


A 

R 


3 


s 


R 


s 


2Rs  +  Rs 
Rg  R 

d(6Vs) 


dV 
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d(6Vn) 

“  dV 

and 

a  =  ( dl/dV) s 

(di/dV)n 

the  ratio  of  dynamic  conductance  in  the  superconducting 
phase  to  the  dynamic  conductance  in  the  normal  phase. 

Then  it  follows  from  (3.7)  that: 

G  =  (l  +  A)  [l  +  ( l-+-pi )  Tj  s  1  [l  +  ( 1+1-0  ( l+^)Tln)  •  (3.8) 


Whenever  the  nature  of  the  I-V  characteristic  of  the  tunnel 
junction  being  measured  permitted,  the  experiment  was 
carried  out  using  R  =  R  .  and  thus  A  =  0.  This  simplifies 

o  II 

the  evaluation  of  G,  since  for  our  specimens  r\s  and  t] 
were  small  (  =  10  J)  while  \i  was  usually  about  10  or  less. 
For  this  special  case  we  have: 


G 


1  + 


( 1+h) 


d(6Vs-6Vn) 

dV 


(3.9) 


In  other  cases  when  Rg  ^  Rr  and  an  uncertainty  of 

AG  =  GAric,  could  be  tolerated  in  G  we  used  the  expression 

G  =  (1  +A)  [1  +  (1+p)  ]  .  (3.10) 


Fortunately  we  could  keep  A  =  0.01  in  these  cases  and 

then  (3.10)  was  accurate  enough.  Where  greater  precision 
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one 

is  required  and  A  ^  0/must  calculate  G  by  (3*8)  which 

is  a  tedious  procedure  indeed. 

In  order  to  show  more  clearly  the  nature  of 

our  bridge  measurements  we  give  the  following  simple 

example.  Figy  3*^  (a)  shows  a  non-linear  current- voltage 

characteristic,  as  well  as  the  voltage,  V* ,  across  R^ . 

The  6Y-Y  curve  appears  in  Fig.  3.^  (b) .  Now  if  we  let 

R  =  R,  »  R-,  and  R  »  R-,  then: 
a  b  d  g  d 


V  =  iR,  -  <5V 
d 

and 

ay  =  R  _  djjvi 

di  d  di 

R  _  d(6V)  OV 
d  dV  di 

Thus 

di  =  1  +  V 

dV  R^ 


which  is  identical  to  (3.7)  under  similar  assumptions. 


C.  Data  Processing 

All  curves  were  traced  with  increasing  voltage 
and  also  with  decreasing  voltage.  This  assured  us  that  there 
was  no  appreciable  phase  shift  in  the  galvanometer 
amplifier  at  the  sweep  rates  we  used  and  any  spurious 


. 
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FIGURE  3.4 


Schematic  showing  how  a  6V-V  curve  (b)  is  obtained 
from  the  current-voltage  curve  (a) . 


(a) 


signals  which  affected  the  two  traces  differently  could 
be  distinguished  from  stable  specimen  characteristics. 

Any  effects  due  to  seismic  disturbances  in  the  neighbour¬ 
hood  of  the  galvanometer  amplifier.,  or  to  heating  in  the 
tunnel  junction,  were  readily  detected  in  this  way.  Only 
those  events  present  in  both  traces  at  the  same  voltage 
were  used  in  the  data  analysis. 

The  curves  were  recorded  on  negative  film* 
using  an  oscilloscope  camera.  The  negatives  were  projected 
through  an  enlarger  onto  a  sheet  of  graph  paper.  Linear 
magnifications  as  high  as  five  were  used.  Two  negatives, 
one  taken  in  the  normal  state  and  one  in  the  superconduct¬ 
ing  state,  having  the  same  scales,  were  traced  on  one  sheet 
of  graph  paper.  The  separation  between  two  such  curves, 
along  the  ordinate  direction,  is  proportional  to  ( 6V  -  6V  ) . 
This  separation  was  measured  as  a  function  of  junction 
voltage  with  as  many  as  500  points  being  obtained  on  some 
of  our  experimental  runs.  The  relative  conductance  was 
then  obtained  easily  when  the  separations  were  multiplied 
by  the  appropriate  scale  factor  to  aUcwfor  the  magnifica¬ 
tion.  This  calculation  was  carried  out  using  the  IBM  1620 
digital  computer.  We  were  able  to  adjust  our  sensitivity 


* 


Kodak  Contrast  Process  Ortho  sheet  film  was  used  and 
developed  in  Kodak  D-8  developer  for  about  45  seconds. 
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for  6V  up  to  a  maximum  of  about  1  jiY / cm  on  the  oscilloscope 
screen,  whereas  the  sensitivity  for  V  was  variable  up  to 
200  jiY / cm  which  was  ample.  In  some  instances  as  many  as 
thirty-five  photographs  were  taken  to  obtain  a  relative 
conductance  curve . 


CHAPTER  IV 


DESIGN,  CONSTRUCTION,  AND  OPERATION  OF  A  He3  CRYOSTAT 


A.  Introduction 

R  4 

The  use  of  He  instead,  of  He  to  obtain  low 
temperatures  by  allowing  the  liquid  helium  to  evaporate 

*)f 

under  its  saturated  vapour  pressure '  has  two  distinct 

advantages.  (i)  As  can  be  seen  in  Fig.  4.1,  the  vapour 

R  4 

pressure  of  He  is  very  much  higher  than  that  of  He  ,  thus 

much  smaller  vacuum  pumps  may  be  used.  (ii)  In  the  temp¬ 
erature  region  of  interest  (above  0.25°K)  He^  does  not 

exhibit  superfluidity,  thus  the  phenomenon  of  film  flow 

4 

which  gives  rise  to  a  large  heat  influx  to  a  He  bath  at 

temperatures  below  2.2°K  is  absent.  Such  heat  influx 

arises  as  a  direct  consequence  of  the  superfluidity  of 
4 

He  .  It  is  due  to  the  rapid  flow  of  a  thin  film  of 
4 

superfluid  He  from  the  bath  along  the  walls  of  the  pumping 
tube  to  regions  of  higher  temperature.  This  film  evapor¬ 
ates  and  has  to  be  pumped  away,  otherwise  it  will  condense 
again  on  the  surface  of  the  bath  where  it  will  deliver 
its  heat  of  vapourization.  These  reasons  make  the  use  of 
He  Impractical  as  a  refrigerant  below  1°K. 


* 


The  pressure  Is  reduced  below  atmospheric  by  the  use 
of  a  vacuum,  pump.  The  temperature  is  then  varied  by 
adjusting  the  pressure  over  the  boiling  liquid. 
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FIGURE  4.1 


Vapour  Pressures  of  the  Stable  Helium  Isotopes. 

Data  for  the  He  scale  was  obtained  from  Sydorlak 

.  .  4 

and  Roberts  (1957).  Data  for  the  He  scale  was 

obtained  from  van  Di jk_,  Durieaux,,  Clement,  and 

Logan  (1958). 
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Temperatures  below  1°K  could  have  been  attained 
by  the  use  of  adiabatic  demagnetization  of  a  paramagnetic 
salt.  However  a  He  cryostat  has  the  advantage  of  sim¬ 
plicity  in  that  no  large  electromagnet  is  required.  For 
experiments  in  which  the  magnetic  state  of  the  specimen 
prior  to  cooling  is  of  importance,,  magnetic  cooling  has 
an  added  disadvantage.  For  example  if  work  is  being 
carried  out  on  superconductors  some  magnetic  flux  is 
trapped  in  the  specimen  during  the  cooling  process.  Thus 
some  type  of  thermal  switch  must  be  incorporated  in  the 
cryostat  in  order  to  isolate  the  specimen  from  the  salt 
pill.  The  specimen  must  then  be  warmed  to  some  temperature 
above  the  superconducting  transition  temperature  to  expel 
the  trapped  flux.,  and  then  allowed  to  cool  again.  For 
these  reasons  the  use  of  a  He  cryostat  has  a  decided 
advantage  over  other  methods  of  refrigeration  when  inves¬ 
tigating  properties  of  superconductors  below  1°K. 

B.  Sources  of  He 

q 

The  natural  abundance  of  He  is  extremely  low, 

_  q  _  5 

constituting  only  10  J  to  10  of  helium  gas  obtained, 
from  natural  gas  wells  ( cf .  Aldrich  and  Nier  1948). 
Commercially  available  He^  is  obtained  artificially  as 
a  by-product  of  tritium  production  in  nuclear  reactors. 


, 
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Li^  +  n 

-*  He4  + 

H3 

4  8 

He  +  HJ 

-  He3  + 

He 

As  the  above  equations  show,  even  with  artificially 

produced  He  ,  one  must  still  separate  the  two  isotopes. 

Early  methods  of  separation  involved  thermal  diffusion. 

Daunt,  Probst  and  Johnston  (1947  and  1948)  showed  that 

He  does  not  participate  in  superfluid  flow  and  thus  showed 

the  way  to  low  temperature  methods  for  the  separation  of 

the  stable  helium  isotopes.  Many  methods  of  separation 

have  since  been  devised  and  are  described  by  Peshkov  and 

Zinov’eva  (1959).  At  the  present  time  the  only  producer 
8 

of  He  on  this  continent  is  the  U.  S.  Atomic  Energy 
Commission.*  The  consumption  of  He^  for  cryogenic  applica¬ 
tions  has  been  greatly  increased  in  the  past  two  years  due 

O 

to  a  reduction  in  the  price  of  HeJ  from  about  $1500  to  $150 
per  liter  (N.T. P. ) . 

Commercial  He  contains  small  traces  of  H  ,  a 
beta  emitter  with  a  half  life  of  12  1/2  years,  and  for 
this  reason  some  precautions  should  be  taken  in  the  initial 
handling  of  He  .  Since  the  beta  particles  are  of  low 
energy  ( 18  keV)  no  shielding  other  than  that  provided  by 
the  gas  container  and  the  metallic  or  glass  pumping  tubes 
of  the  cryostat  need  be  used.  However  precautions  should 
be  taken  to  prevent  the  possible  ingestion  of  tritium, 

*  He3  is  distributed  through  contract  with  the  Monsanto 
chemical  Company,  Mound  Laboratory,  Miamisburg,  Ohio. 
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e.g.  glassblowing  on  sections  of  a  cryostat  which  have 
been  in  contact  with  H  should  be  avoided. 

Q 

C.  General  Features  of  He  Cryostats 

A  number  of  He  cryostats*  have  been  described 

in  the  literature.  In  all  cases  the  principle  of  their 

operation  is  very  similar.,  consisting  of  a  He  stage  that 

4 

operates  within  a  He  cryostat  which  has  been  cooled  to 
^  1.25°K.  The  He^  is  liquefied  in  a  condenser  at  1.25°K 
and  the  liquid  is  then  allowed  to  enter  a  chamber  (He^ 
chamber) .  The  experimental  specimen  is  placed  in  thermal 
contact  with  the  He  chamber.  Variation  of  the  pressure 
above  the  liquid  He  allows  one  to  control  the  temperature 
of  the  specimen.  Two  types  of  these  cryostats  have  been 
built.  (i)  The  "single  shot"  cryostats  in  which  the  He^ 
is  allowed  to  condense  in  the  He^  pumping  tube  until 
sufficient  liquid  has  accumulated  in  the  He  chamber. 

The  temperature  may  then  be  reduced  by  pumping  on  the 
He  .  After  all  the  He  is  pumped  away  it  is  necessary  to 
recondense  it  again  at  1.25°K.  (ii)  The  "refrigerator" 
in  which  the  He  is  condensed  in  a  special  condensing 
capillary  and  is  allowed  to  enter  the  He^  chamber  con- 


Roberts  and  Sydoriak  (1955)*  Seidel  and  Keesom  (1958)* 
Reich  and  Garwin  (1958),  Ambler  and  Dove  (1961),  Martin 
(1961),  de  Bruyn  Ouboter,  Taconis,  LePair,  and  Beenakker 
(i960),  van  den  Meijdenberg,  Taconis,  and  Le  Pair  (1961), 
and  Zinov'eva  and  Peshkov  (1959)- 
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3  Q 

tinuously,  i.e.  while  the  He^  is  being  pumped.  The  HeJ 

which  is  pumped  away  is  recondensed  in  a  cyclic  manner. 

This  type  of  cryostat  allows  maintenance  of  a  particular 

temperature  for  prolonged  periods.  Because  it  is  cyclic 

it  requires  smaller  amounts  of  He  ,  since  the  coolant  is 

continuously  replenished.  Due  to  the  price  and  scarcity 

of  He  it  is  always  recovered  and  re-used.  This  puts  a 

very  stringent  requirement  on  the  He  system  in  that  it 

has  to  be  free  even  of  " superleaks" ,  which  are  leaks  that 

are  too  small  to  allow  ordinary  fluids  to  pass,,  but  permit 
4 

superfluid  He  to  pass.  Such  leaks  would  contaminate  the 

7  7  4 

He  and  if  they  were  large  enough  the  He  -He  mixture 

would  exhibit  a  lambda  transition,  the  temperature  of 

4 

which  increases  with  He  concentration.  Such  mixtures 

tend  to  limit  the  ultimate  temperature  attainable  with 
■=>  ^ 

a  He  cryostat.  Our  He  cryostat  was  designed  in  such  a 
way  that  it  served  both  as  a  refrigerator  and  a  single 
shot  cryostat. 

D.  Description  of  the  He  Cryostat 

A  complete  schematic  diagram  of  the  He  cryostat 
and  auxiliary  systems  is  given  in  Fig.  4.2.  A  set  of 
machine  drawings  of  the  cryostat  appears  in  Appendix  III. 
The  experimental  chamber  of  the  cryostat  is  enclosed  by 
a  brass  can  K.  This  can  is  soldered  to  a  brass  plate,  D, 


FIGURE  4.2 


Schematic  diagram  of  the  He 


3 


Cryostat . 


manometer  reference  vacuum  line 
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using  Wood's  alloy *  so  that  it  can  readily  be  removed 

for  changing  specimens.  The  experimental  chamber  can  be 

evacuated  through  two  pumping  tubes,,  L.  These  tubes  are 

made  of  3/8  inch  diameter  stainless  steel  having  a  wall 

thickness  of  0.004  inch.  The  reason  for  using  such  thin 

4 

walled  tubing  is  to  minimize  the  heat  influx  into  the  He 

in  the  dewar .  The  vacuum  inside  the  experimental  chamber 

serves  to  thermally  isolate  its  contents  from  the  surround- 
4 

ing  He'  dewar  (the  helium  in  this  dewar  is  allowed  to  boil 

at  atmospheric  pressure,  so  that  T  =  4°K) .  The  electrical 

leads  are  brought  to  the  experimental  chamber  through 

one  of  these  pumping  tubes.  They  enter  the  pumping  tube 

through  a  Kovar  to  glass  seal  at  the  top  of  the  cryostat. 

4  i  o 

These  wires  are  thermally  anchored  to  the  He  bath  at  4  K. 

This  anchoring  is  achieved  by  tightly  wrapping  several 

feet  of  wire  around  a  copper  post  which  is  soldered  to  D. 

4 

The  wires  are  also  thermally  anchored  to  the  He  chamber, 

F,  with  an  intervening  length  of  about  one  foot  of  wire 
freely  suspended  in  the  vacuum  jacket.  Lengths  of  manganin 
wire  which  have  been  tinned  with  50/50  tin-lead  solder 
linked  these  wires  to  the  specimen  or  carbon  resistance 
thermometers.  This  provides  a  superconducting  link  of 
low  thermal  conductivity.  Small  pieces  of  blackened 
copper  foil,  C,  serve  as  radiation  traps  for  the  two 
pumping  tubes.  Provision  is  made  on  the  vacuum  system 
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to  admit  low  pressure  He  gas  to  the  experimental  chamber. 

This  exchange  gas  is  used  to  achieve  thermal  equilibrium 

4 

between  the  liquid  He  in  the  dewar  and  the  contents  of  the 
specimen  chamber. 

4 

Chamber  F  contains  liquid  He  .  The  chamber  is 
filled  through  a  combination  needle  valve  and  radiation 
trap,,  A,  a  detailed  view  of  which  appears  in  Fig.  A-2, 
Appendix  III.  A  short  length  of  Tygon  tubing  attached  to 
A  and  extending  several  inches  below  plate  D  allows  one 
to  fill  chamber  F  even  if  the  liquid  helium  level  is  below 
plate  D.  The  temperature  of  chamber  F  can  be  reduced  from 
4°K  to  1.25°K  by  pumping  with  the  He^  pump.  Intermediate 
temperatures  can  be  maintained  by  using  a  Cartesian  mano- 
stat  of  the  type  described  by  Gilmont  (1951).  The  He^ 
which  has  been  pumped  from  chamber  F  passes  from  the 
exhaust  of  the  pump  through  cotton  wool  filter  2  into 
the  helium  recovery  system.  An  estimated  25  cubic  feet 

of  He^  gas  (N.T.P.)  are  recovered  from  this  chamber  each 

4  ,n 

run.  The  He  chamber  is  connected  by  means  of  a  l/o 

inch  Inconel  tube  to  a  mercury  and  an  oil  manometer. 

These  manometers  are  used  to  determine  the  vapour  pres- 

4 

sure  from  which  the  temperature  of  the  He  bath  can  be 
deduced.  The  thermomolecular  pressure  correction  had  to 
be  applied  only  at  the  lowest  temperatures  (T  =  1.25°K) 


' 
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and  was  of  the  order  of  1%.  For  this  correction  the 

values  of  Roberts  and  Sydoriak  (1956)  were  used. 

v 

The  He^  part  of  the  system  was  built  around  a 
Welch  model  14-02  KBG  vacuum  pump.  This  is  a  rotary  pump 
with  a  speed  of  about  70  liters  per  minute  at  1  p  pressure. 
It  has  a  specially  designed  shaft  seal  ( cf .  Reich  and 
Garwin  1958)  whose  leak  rate  is  about  200  times  lower 

than  that  of  a  conventional  pump.  Such  a  modified  pump 

3 

is  necessary  to  avoid  contamination  or  loss  of  the  He  . 
o 

The  He^  (1.12  liters  N.T.P.)  is  stored  partly  in  the  pump 
body  and  partly  in  a  brass  storage  vessel  at  a  pressure 
of  about  2/3  of  an  atmosphere.  Low  pressure  storage  is 
used  in  order  to  minimize  potential  loss  due  to  small 
leaks.  Any  air  collected  through  such  small  leaks  can 
readily  be  separated  from  the  He  by  freezing  and  then 
released  later.  The  exhaust  side  of  the  He  pump  plays 
the  role  of  the  compressor  in  the  refrigerator  system. 

The  compressor  pressure  is  always  below  one  atmosphere, 
and  about  equal  to  the  saturated  vapour  pressure  of  the 
He^  at  the  condensation  temperature  (T  =  1.25°K).  Filter 
1  (glass  wool)  collects  oil  vapours  from  the  He'  gas  as 
it  leaves  the  pump.  Condensation  takes  place  by  letting 
the  HeJ  gas  pass  through  valves  2  and  3  into  a  1.0  milli¬ 
meter  capillary  (condensing  line).  After  passing  through 
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i  o 

the  He  bath  at  4  K  the  capillary  narrows  to  a  diameter 
of  0.5  millimeters  and  feeds  coil  E  which  consists  of  a 
one  meter  length  of  0.5  millimeter  capillary  wound  into  a 

coil  with  a  diameter  of  1/8  inch.  This  coil  is  located  in 

4  o 

the  He  chamber  and  is  at  a  temperature  of  1.25  K.  The 

3 

He  liquifies  inside  this  coil  provided  the  compressor 
pressure  is  25  millimeters  of  mercury  or  more  (vapour 
pressure  of  He^  at  1.25°K).  The  He^  thus  liquified  passes 

3 

through  a  porous  plug.,  G-,  into  the  He  chamber.  This 

3 

plug  serves  to  limit  the  injection  rate  of  He  into  the 
chamber.  Too  rapid  an  injection  rate  will  raise  the 
ultimate  temperature  because  flow  resistance  of  the  tubes 

3 

connecting  the  He  chamber  to  the  pump  finally  limits  the 
refrigeration  rate.  The  plug  is  made  by  inserting  a 
piece  of  porous  stainless  steel  which  has  been  ground 
down  to  0.050  inch  diameter  into  a  copper  tube  having  an 

/  r  * 

inside  diameter  of  1/16  inch.  The  copper  tube  was  then 

drawn  down  until  it  fitted  tightly  around  the  stainless 

steel.  The  plug  used  in  this  cryostat  had  a  flow  rate 
4  3 

for  He  gas  of  about  1.6  cm  /min  under  a  pressure  differ¬ 
ence  of  1/2  atmosphere  at  room  temperature.  The  specimen 

3 

is  attached  to  post  J,  which  together  with  the  He  chamber 
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is  machined  from  one  piece  of  copper  (see  Appendix  III). 

Q 

Some  control  of  the  temperature  of  the  He  chamber  is 

obtained  by  reducing  the  pumping  speed  with  valve  4.  In 

this  way  the  cryostat  operates  in  a  continuous  cycle  as  a 

refrigerator.  Single  shot  operation  is  achieved  by  closing 

valves  2  and  3  and  opening  valve  1  so  that  the  He^  condenses 

down  the  pumping  tube.  This  temporarily  raises  the  temper- 

4 

ature  of  the  He  chamber  but  when  equilibrium  is  re-established 

o 

valve  1  is  closed  and  valve  4  opened  to  cool  the  He  .  A 
system  for  leak  detection  preceding  every  run  is  also  incor¬ 
porated  in  the  apparatus.  The  details  of  this  system  are 
discussed  in  Section  G-. 

E.  Discussion  of  Some  of  the  .Design  Parameters  of  the 

q 

Liquid  He  Cryostat 

Due  to  the  scarcity  and  price  of  He  at  the  time 
the  cryostat  was  designed,  we  adopted  the  basic  philosophy 
of  keeping  the  number  of  joints,  particularly  at  the  low 
temperature  end  of  the  cryostat,  to  a  minimum.  This  has 
resulted  in  a  performance  which  is  somewhat  less  than 
optimum,  in  particular  the  ultimate  attainable  temperature 
is  only  0.47°K.  Modifications  to  remedy  this  situation 
are  proposed  in  Section  H.  In  the  present  section  we 
concern  ourselves  only  with  the  design  parameters  of  the 
existing  system.  The  design  problems  encountered  in 
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3 

building  a  He-^  cryostat  are  in  many  respects  similar  to 
those  of  constructing  an  adiabatic  demagnetization  cryostat. 
The  major  fundamental  difference  between  the  two  systems 
is  that  the  low  temperature  part  of  the  apparatus  (salt 
pill)  in  a  demagnetization  cryostat  has  a  much  higher  heat 
capacity,,  than  does  a  He  chamber  and  its  content.  This 
results  in  He  systems  having  a  much  more  rapid  equilibrium 
time  than  that  of  demagnetization  systems,,  as  well  as 
greater  temperature  instabilities  than  demagnetization 
cryostats  under  similar  heat  influxes. 

o 

(i)  Pumping  Speed  of  the  HeJ  system: 

The  equations  governing  the  flow  of  gases  through 
pipes  are  discussed  in  detail  by  Dushman  (1949).  There  are 
however  a  number  of  problems  peculiar  to  cryogenic  equip¬ 
ment  which  make  difficult  accurate  estimation  of  the  size 
of  pumping  tubes  required  for  particular  applications.  In 
an  isothermal  situation  the  equations  of  gas  conductance 
can  be  solved  exactly  in  two  limiting  cases.  The  first 
case.,  free  molecular  flow.,  occurs  when  the  mean  free  path 
of  the  molecules  is  larger  than  the  lateral  dimensions  of 
the  gas  passage.  In  this  case  conduction  is  determined 
entirely  by  the  geometry  of  the  tubes.  The  second  case, 
Poiseuille  flow,  arises  when  the  mean  free  path  of  the 
molecules  is  much  less  than  the  lateral  dimensions  of  the 
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gas  passage.  The  pumping  speed  then  depends  on  the  gas 
viscosity,  pressure  gradient,  and  geometry  of  the  pumping 
tube.  Even  in  these  two  cases  exact  calculations  can 
only  be  made  for  pumping  tubes  whose  lengths  are  about 
fifty  times  their  diameter  (infinite  tubes)  or  for  orifices 
(zero  length).  In  practical  applications  of  the  type 
considered  here  neither  of  these  cases  is  applicable  and 
the  corrections  which  must  be  made  for  end  effects  are 
large.  Over  the  pressure  and  temperature  ranges  involved 
the  type  of  flow  is  neither  viscous  nor  free  molecular 
flow,  and  unfortunately  no  satisfactory  interpolation 
formula  between  these  two  types  of  flow  exists.  An 
added  difficulty  arises  due  to  the  fact  that  most  of  the 
pumping  tubes  are  under  a  temperature  gradient.  In  cal¬ 
culating  the  required  pumping  tube  diameter  on  the  basis 
of  free  molecular  flow  for  an  ultimate  temperature  of 
0.35°K,  a  diameter  of  about  1/4  inch  is  sufficient  for  the 
section  of  the  He^  pumping  tube  between  the  4°K  bath  and 
the  HeJ  chamber.  It  is  further  found  that  the  2  1/2  inch 

diameter  section  of  the  pumping  tube  at  room  temperature 

q 

limits  the  ultimate  pumping  speed  of  the  whole  He^  pump- 
ing  line.  The  section  of  the  He^  pumping  line  between 
4°K  and  the  He^  chamber  was  made  3/4  inch  in  diameter  in 
order  to  avoid  extra  solder  joints.  This  results  in  a 
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large  heat  leak  to  the  He  bath  and  limits  the  ultimate 
temperature.  Because  of  this  it  is  impossible  to  maintain 
intermediate  temperatures  (between  0.47°K  and  1.25°K) 
using  a  single  shot  process.  Remedies  for  this  situation 
are  given  in  section  H. 

(ii)  The  dead  volume: 

In  order  to  make  the  most  efficient  use  of  the 
available  He  ,  the  dead  volume,  that  is  the  volume  which 
contains  He  gas  when  the  maximum  amount  has  been  condensed, 
must  be  kept  to  a  minimum.  For  this  reason  the  volume  of 
the  condenser  side  of  the  cryostat  was  made  of  capillary 
tubing.  Most  of  the  dead  volume  is  in  the  pumping  tubes 
and  this  is  difficult  to  reduce  without  impairing  the 
pumping  speed.  The  distribution  of  He  by  weight,  when 
all  the  condensable  He^  is  in  liquid  form  at  1.25°K  is  shown 


below: 

Amount  at  room  temperature  6$ 

Amount  between  295°K  and  4°K  2$ 

Amount  at  4°K  8$ 

Amount  at  1.25°K  but  not  in 

the  He^  chamber  18 $ 

Amount  in  He^  chamber  66$ 
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(iii)  Heat  Influx: 

Estimates  of  the  heat  flow  to  the  He  chamber 

were  made  assuming  the  temperature  to  be  below  0.5°K  and 
4  n 

the  He  bath  to  be  at  1.25  K.  The  biggest  source  of  heat 

influx  to  the  HeJ  arises  from  the  conduction  by  He  gas 
4  3 

between  the  He  and  He  baths.  This  mechanism  can  contri¬ 
bute  as  much  as  40  microwatts;  this  upper  limit  is 
obtained  by  calculating  the  conduction  of  a  static  gas 

O 

column.  The  3/4  inch  pumping  tube  joining  the  He3  and 
4 

He  baths  and  the  electrical  leads  conduct  about  15 
microwatts  each.  Initially  the  heat  leak  through  the 
electrical  leads  was  very  large  due  to  improper  anchoring 
of  the  wires  at  4°K  and  at  1.25°K.  Estimates  of  the 
heat  leak  due  to  inefficient  anchoring  obtained  from 
observations  of  the  warming  rate  of  the  cryostat  show  this 
contribution  to  be  about  2  microwatts. 

A  major  source  of  heat  influx  is  mechanical 
vibration.  This  is  minimized  to  some  extent  by  mounting 
the  cryostat  on  a  pillar  which  is  mechanically  de-coupled 

3 

from  the  floor.  The  He  pump  Is  connected  to  the  cryostat 

4 

by  means  of  a  bellows,  while  the  He  pump  is  de-coupled 
using  flexible  vacuum  hose.  Despite  such  precautions 
slight  mechanical  vibrations  are  transmitted  to  the 
cryostat  and  result  in  additional  heat  influx. 


, 
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(iv)  The  Kapitza  Resistance: 

Kapitza  (19^1)  showed  that  in  the  case  of 
4  o 

liquid  He  below  about  1  K  there  is  a  temperature  dis¬ 
continuity  at  the  boundary  between  the  liquid  helium  and 
the  solid  surface  of  the  container.  At  the  present  time 
no  proper  theory  of  this  effect  exists.  The  experiments 
of  Fairbank  and  Lee  (1957)  show  that  this  phenomenon  also 
occurs  with  He  .  Their  results  indicate  that  with  He 
the  thermal  resistivity,,  W,  is  given  by 

¥  =  130/T2  cm2  deg  K  watt  \ 

where  T  is  the  absolute  temperature  at  the  boundary.  The 
temperature  discontinuity.,  AT,  is  then  given  in  terms  of 
the  rate  of  heat  transfer.,  in  watts  across  the  boundary 
by 

AT  =  130  Q/AT2  , 

2 

where  A  is  the  surface  area  in  cm  across  which  heat 

q 

transfer  takes  place.  The  He  chamber  was  designed 
to  have  a  surface  area  which  is  large  enough  that  this 
effect  is  not  too  important.  For  a  heat  load  of  about 
1C>  watts  the  maximum  temperature  discontinuity  is 
calculated  to  be  three  millidegrees  K. 


' 
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F.  Thermometer  Calibration 

Two  types  of  commercial  carbon  resistors  were 
used  as  resistance  thermometers,,  one  was  a  10  ohm  Allen- 
Bradley  resistor,,  the  other  a  220  ohm  ,5peer.  The  resis¬ 
tances  of  these  thermometers  were  determined  using  an 
isolating  potential  comparator  similar  in  principle  to 

that  described  by  Dauphinee  (1953).  These  thermometers 

4 

were  calibrated  against  the  He  vapour  pressure  scale  at 
several  temperatures  between  4.2°K  and  1.25°KJ(  and  an 
additional  comparison  was  made  with  the  He^  vapour  pres¬ 
sure  scale  at  1.25°K.  The  range  from  1 . 25°K  to  0.47°K 
was  covered  using  an  extrapolation  of  the  three  constant 
equation  of  Markham,,  Netzel  and  Dillinger  \  1957 )  *• 

InR/T  =  a(lnR)2  +  b(lnR)  +  c 

Here  R  is  the  thermometer  resistance,  T  is  the  absolute 
temperature,  and  the  constants  a,  b,  and  c  are  determined 
by  a  least  squares  fit  to  the  calibration  points.  In 
order  to  determine  the  validity  of  such  an  extrapolation 
two  different  thermometers  were  attached  to  the  specimen 
post  and  calibrated  in  the  manner  described  above.  At 
the  same  time  a  piece  of  high  purity  cadmium  wire  was 
also  attached  to  the  specimen  post.  The  cadmium  had  both 
current  and  potential  leads  attached  to  it.  The  resistance 


* 
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of  the  cadmium  was  determined  as  a  function  of  temperature. 
The  resistance  measurements  were  made  by  comparing  the 
voltage  drop  across  a  standard  resistor  in  series  with  the 
specimen  to  the  voltage  drop  across  the  specimen.  The 
voltage  measuring  system  was  that  described  by  Adler  and 
Woods  (1962).  The  temperature  at  which  the  cadmium  became 
superconducting  was  determined  with  both  resistance  thermo¬ 
meters.  The  superconducting  transition  temperature  of 
cadmium  obtained  by  extrapolation  of  the  two  thermometers 
was  0.52  +  0.015°K.  This  is  in  good  agreement  with  the 
value  of  0.52  -  0.01°K  obtained  calorimetrically  by 
Martin  (1961a)  for  cadmium  of  similar  purity.  The  reason 
that  the  temperature  scale  was  not  established  using  the 
He  vapour  pressure  scale  is  that  this  would  have  required 
at  least  one  additional  liter  of  He^  (N.T.P.)  for  filling 
a  separate  manometer  system. 

CL  Operating  Procedure  of  the  He^  Cryostat 

This  type  of  cryostat  differs  considerably  from 
4 

conventional  He  cryostats  in  that  a  rather  labourious 
leak  testing  procedure  must  be  carried  out  in  the  process 
of  cooling  down  the  apparatus.  Leak  testing  is  performed 
using  a  Veeco  leak  detector.  This  instrument  is  a  mass 

4 

spectrometer  which  is  sensitive  only  to  He  ,  and  is  capable 
of  detecting  very  small  amounts  of  this  gas  (2  x  10  ^std 
cm^/sec)  leaking  into  the  system  to  which  it  is  attached. 


. 
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Before  every  run  the  He  chamber,  the  experimental 

3  3 

chamber,  and  the  He  system  are  evacuated.  The  HeJ  system 

is  evacuated  through  the  auxiliary  vacuum  system  shown  in 

Pig.  4.2.  This  auxiliary  system  is  used  in  order  to  avoid 

contaminating  the  leak  detector's  vacuum  system.  When  the 

1  4 

pressure  in  the  HeJ  system  has  been  reduced  to  below  10  mm 

of  Hg  (black-out  on  discharge  gauge,  M)  valve  7  is  opened. 

4 

He  from  a  cylinder  is  sprayed  on  all  the  exposed  parts  of 

3 

the  He  system.  If  the  leak  indication  then  remains  below 
-9  3  / 

about  10  std  cm  /sec  the  system  is  assumed  to  be  leak 
4 

tight.  The  He  chamber  is  then  filled  with  about  1/2 

4  4 

atmosphere  of  He  gas  from  a  cylinder,  and  the  He  system 

is  sealed.  If  no  increase  in  the  leak-rate  is  evident 

3 

after  about  five  minutes  the  He  condensing  coil  and  length 

of  He^  pumping  tube  are  assumed  to  be  leak  tight.  A  similar 

procedure  is  carried  out  filling  the  experimental  chamber 

with  He^  exchange  gas  to  a  pressure  of  about  10-^  mm  of  Hg . 

This  tests  the  He^  chamber  and  porous  plug  for  leaks. 

4 

With  the  He  in  the  two  chambers  as  described  above  the 
apparatus  is  then  cooled  down  to  liquid  air  temperatures. 

If  no  leak  indication  appears  in  the  ensuing  hour  the 
apparatus  is  then  cooled  down  to  liquid  helium  temperature. 
The  He^  in  the  dewar  serves  to  test  the  pumping  tubes  for 


leaks.  When  the  liquid  helium  level  in  the  dewar  is 


. 


j  a 


sufficiently  high  needle  valve.  A,  is  opened  several 

4 

turns.  This  admits  liquid  helium  into  the  He  chamber. 

4  / 

When  the  pressure  in  the  He  chamber  (as  monitored 

either  by  the  manometers  or  by  the  Bourdon  gauge)  reaches 

4 

atmospheric  pressure,  the  vent  near  the  He  pump  may  be 

opened.  At  this  point  if  no  increase  in  the  leak  rate 

is  noted  the  experimental  chamber  is  evacuated  to  a 

—  (5 

gauge  pressure  of  about  10  mm  of  Hg.  This  takes  about 

one  hour  to  one  and  one  half  hours .  If  the  leak  rate 

4 

has  remained  unaltered,  the  vent  near  the  He  pump  and 
the  needle  valve  may  be  closed  and  the  system  pumped  down 
to  1.25°K.  At  this  temperature  superleaks  in  the  He^ 
system  should  become  evident.  If  none  appear  the  leak 
detector  is  removed  from  the  system. 

Before  condensing  the  He  through  the  capillary 
it  is  allowed  to  condense  down  the  pumping  tube  by  open¬ 
ing  valve  1.  A  few  minutes  later  it  is  pumped  off  again. 
This  process  must  be  carried  out  several  times  at  the 

3 

beginning  of  every  run.  The  reason  for  this  is  that  He 
system  collects  small  amounts  of  air  between  runs.  This 
condensation  process  freezes  the  air  in  the  large  pumping 
tube  where  it  can  not  create  a  blockage.  The  air  collected 
in  the  large  pumping  tube  may  be  released  at  the  end  of  the 
run  through  the  He^  manometer  system.  During  cyclic  oper¬ 
ation  the  temperature  is  regulated  by  adjusting  the  ten 
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turn  metering  valve  (valve  3)  in  the  condensing  line,  and 

also  adjusting  the  pumping  speed  using  valve  4.  During 

"2 

single  shot  operation  the  He  is  condensed  through  the 

large  pumping  tube,  and  a  final  temperature  of  0.47°K  may 

be  obtained  with  valve  4  wide  open.  It  is  not  possible 

to  maintain  intermediate  temperatures  in  the  single  shot 

process  because  the  heat  influx  down  the  3/4  inch  pumping 

tube  leading  to  the  He  chamber  is  too  large  at  all  but 

the  lowest  pressure.  Temperatures  between  4°K  and  1.25°K 

4 

are  obtained  by  regulating  the  pressure  of  the  He  chamber 

using  the  cartesian  manostat;  the  specimen  is  then  linked 

4  3  3 

thermally  to  the  He  bath  by  means  of  the  He  in  the  He 

pumping  tube. 

* 

H.  Proposed  Modifications  to  the  He^  Cryostat 

At  the  time  the  He^  cryostat  was  originally 

designed,  allowances  were  made  so  that  future  modifications 

might  be  incorporated  without  undue  difficulty.  A  schematic 

diagram  of  a  proposed  modification  of  the  cryostat  is  shown 

in  Pig.  4.3.  This  modification  would  reduce  the  influx  to 

the  HeJ  bath,  reduce  the  dead  volume  by  l6  percent,  and 

provide  for  good  heat  exchange  between  the  incoming  and 

evaporating  He^.  The  latter  is  achieved  by  placing  the 

2 

condensing  coil  in  the  He  pumping  tube.  In  the  modified 

3 

cryostat  the  condensing  capillary  would  enter  the  He 


. 
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FIGURE  4.3 


Schematic  diagram  showing  proposed  modifications 
3 

to  the  He  cryostat.  Although  not  a  scale  drawing 


this  figure  is  approximately  true  size. 


specimen  post 
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•=> 

pumping  line  through  a  ferrule  in  the  He  radiation  trap 
(Pig.  A-3,  Appendix  III).  Both  the  modified  He^  and  the 
He^  chambers  are  made  of  copper.  The  diaphragn  in  the  He^ 
pumping  line  is  similar  to  the  one  in  the  present  cryostat 
and  serves  to  limit  film  flow.  The  pumping  tube  joining 
the  HeJ  and  He  chambers  is  to  be  3/l6  inch  thin  wall  inconel 
tubing.  It  is  expected  that  this  modification  will  enable 
temperatures  of  0.35°K  to  be  reached.,  and  will  allow  the 
cryostat  to  handle  larger  heat  loads. 
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CHAPTER  V 


EXPERIMENTAL  METHOD 

A.  Preparation  of  Tunnel  Junctions 


All  tunnel  junctions  used  in  these  experiments 
involved  thin  evaporated  metal  films.  These  films  were 
deposited  on  1/4"  x  1"  thin  soda  glass  slides.  The  vari¬ 
ous  stages  of  specimen  preparation  are  shown  schematically 
in  Fig.  5.1.  The  glass  slide  (a  portion  of  a  standard 
microscope  slide)  was  thoroughly  cleaned  using  household 
detergent  and  water.  After  drip  drying  it  was  heated  in 
a  Bunsen  flame  to  a  temperature  at  which  the  edges  just 
became  fire  polished.  This  anneals  the  slide  as  well 
as  removing  minor  surface  irregularities.  Four 
copper  wires  were  then  attached  to  the  clean  glass  sur¬ 
face  using  pure  indium  solder  as  shown  in  Fig.  5.1  (a). 

The  copper  wires  were  tinned  with  50/50  tin-lead  solder. 
This  tinning  insured  that  the  specimen  leads  were  super¬ 
conducting  at  liquid  helium  temperatures,,  thus  reducing 
the  heat  input  to  the  tunnel  junction.  The  glass  slide 
was  then  mounted  in  the  evaporator.  The  first  metallic 
layer  (usually  aluminum)  was  then  evaporated  through  an 
aluminum  mask;  Fig.  5.1  (b)  shows  such  a  layer  on  the 


Metal  films  deposited  on  slides  which  had  not  been 
flamed  often  lacked  continuity. 
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FIGURE  5.1 


Schematic  showing  the  various  stages  of  specimen 
preparation . 


(a)  (b) 


(c) 


(d) 
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glass.  Different  kinds  of  filaments  were  tried,,  the  gener¬ 
ally  used  types  being  a  tungsten  helix*  and  molybdenum 

wire  wound  in  a  spiral  to  form  a  basket.  Pressures  of 

_c, 

about  10  'mm  of  Hg  (sometimes  lower)  were  maintained 
during  evaporation.  Care  was  taken  to  let  some  of  the 
material  evaporate  onto  a  shutter  plate  before  any  was 
allowed  to  reach  the  glass  plate.  This  decreased  the 
probability  of  contamination  of  the  metal  film  on  the  glass. 
Evaporation  was  allowed  to  proceed  until  the  filament  was 
no  longer  visible  through  the  glass.  The  resulting  films 
all  had  mirrorlike  finishes.  Air  was  then  admitted  into 
the  evaporation  chamber  to  allow  the  specimen  to  oxidize . 
Oxidation  procedures  varied*  due  to  the  varying  moisture 
content  of  the  air.  When  the  moisture  content  of  the  air 
was  low*  the  air  was  allowed  to  pass  over  water  before 
entering  the  evaporator.  The  specimen  with  both  the  metal 
and  metallic  oxide  layer  is  shown  in  Fig.  5.1  (c).  Finally 
a  second  metallic  layer  was  evaporated*  completing  the 
specimen  as  shown  in  Fig.  5.1  (d). 

If  the  specimen  showed  tunneling  characteristics* 
it  was  mounted  and  used  in  the  experiments.  The  criteria 
used  to  determine  whether  the  major  contribution  to  the 
current  was  due  to  tunneling  was  that  given  by  Fisher  and 
Oiaever  (1961).  The  ratio  of  useless  to  useful  specimens 
was  about  6/1.  Typical  current-voltage  characteristics 


■ 


for  an  A1  -  Al^O^  -  A1  junction  are  shown  in  Pigs.  5.2, 

5*3;  and.  5.5.  Other  tunnel  junctions  behaved  much 

the  same  way  except  that  those  of  type  A  -  AxOy  -  B  (A 
different  from  B)  had  their  center  of  symmetry  displaced 
by  a  few  tenths  of  millivolts  from  the  origin  along  the 
abscissa.  Fig.  5-3  shows  a  comparison  of  I-V  character¬ 
istics  between  our  specimens  and  those  reported  by  Fisher 
and  G-iaever  (1961). 

The  following  types  of  tunnel  junctions  were 

successfully  made  and  used  in  these  experiments : 

Pb  -  AlyO^  -  Al,  Zn  -  Alo0^  -  Al,  A1  -  Al^O^-  Al,  and 

In  -  AlyO^-  Al .  Attempts  were  made  at  using  various 

other  oxide  layers  but  all  of  these  failed.  Among  those 

tried  were:  Cu  -  Cu  0  -  In.  Fe  -  Fe  0  -  In,  In  -  In^O  -  In, 

x  y  ’  x  y  x  y 

Sn  -  AI2O3-  Al,  Pb  -  PbO  -  Pb,  Al  -  Al^O^-  Tl,  and 
Mg  -  MgO  -  In.  One  of  the  difficulties  with  the  magnesium 
was  that  even  though  evaporated  through  a  mask  the  magnesium 
spread  out  and  covered  an  area  much  larger  than  that  of 
the  mask,  and  in  many  cases  the  whole  glass  slide.  This 
occured  even  when  the  mask  was  only  separated  by  about 
2  mm  from  the  glass  slide.  A  possible  cause  for  this 
behaviour  is  the  rather  high  vapour  pressure  of  magnesium 
at  the  temperatures  at  which  we  attempted  evaporation. 

This  is  unfortunate  bacause  magnesium  might  have  provided 
a  metal  substrate  which  does  not  become  superconducting  and 
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FIGURE  5.2 


Current-voltage  curves  for  a  typical  Al-Al.^O^-Al 
tunnel  junction  at  three  different  temperatures. 
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FIGURE  5.3 


Semi-log  plot  of  the  I-V  curves  of  typical 
AI-AI^Cr-AI  tunnel  junctions. 
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FIGURE  5.4 


Linearity  of  a  typical  Al-Al^Cn-Al  tunnel  junction 
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FIGURE  5.5 


Symmetry  of  a  typical  Al-AlpO^-Al  tunnel  junction. 
(Forward,  voltage  is  with  the  aluminum  film  nearest 


the  glass  made  positive) . 
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at  the  same  time  is  similar  to  aluminum  in  its  oxide 
properties.  The  advantage  of  a  normal  metal  substrate 

is  that  measurements  can  then  be  carried  out  with  only 

\ 

one  component  of  the  junction  superconducting  and  then  the 
relative  conductance  gives  the  ratio  of  the  densities  of 
states  of  the  superconducting  to  normal  phases  directly. 

Some  attempts  were  also  made  to  produce  tunnel 
junctions  with  mercury  as  one  of  the  metals.  One  such 
junction  was  made  using  a  "Kovar/  tube  that  had  an  oxidized 
surface,  and  was  subsequently  filled  with  mercury.  The 
oxide  layer  was  too  thick  and  although  conduction  through 
the  "Kovar"  oxide  was  partly  due  to  tunneling,,  the  bulk 
of  the  conduction  was  due  to  a  semiconductor  effect.  This 
mercury- "Kovar"  oxide-"Kovar"  junction  had  a  room  temper¬ 
ature  resistance  of  1/4  ohm,  while  its  resistance  at  1.25°K 
was  about  1/2  megohm.  All  that  this  experiment  demonstrated 
was  that  "Kovar"  oxide  is  a  semiconductor.  Other  attempts 
involved  mechanically  polished  copper,  coated  with  a  very 
thin  layer  of  polyvinyl  formvar  which  floated  in  a  pool  of 
mercury  that  was  subsequently  frozen.  This  type  of  junction 
could  be  made  to  exhibit  tunneling  properties  down  to  79°K, 
but  would  breakdown  after  about  an  hour.  The  reason  for 
the  breakdowns  were  not  ascertained.  Further  attempts  at 
producing  tunnel  junctions  with  mercury  should  be  made 
because  mercury  is  an  interesting  superconductor,  since 
the  electron-phonon  interaction  in  this  metal  is  very  strong. 
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6.  Specimen  Mounting 

The  specimens  on  glass  slides  were  mounted,  one 
at  a  time  in  a  "Kovar"  to  glass  seal  shown  in  Pig.  5-6. 

The  top  cap  of  the  seal  holds  the  specimen,  while  the  cur¬ 
rent  and  potential  leads  are  brought  out  through  four 
small  "Kovar"  tubes  at  the  bottom.  These  small  tubes  are 
sealed  by  allowing  solder  to  flow  into  them.  The  top  cap 
contains  an  indium  "0"  ring  seal  (Vacuum  seal).  There  is 

a  small  tube  in  the  top  cap  which  allows  the  assembly  to 

4 

be  evacuated  and  then  filled  with  He  gas  to  atmospheric 

pressure.  A  few  pellets  of  solder  which  had  been  placed 

in  the  small  pumping  tube  before  evacuation  are  melted 

4 

to  seal  the  assembly.  This  holder  with  the  He  gas  and 

specimen  inside,  was  then  soldered  onto  the  specimen  post 

of  the  cryostat.  The  copper  (about  0.040  inch  diameter) 

joining  the  top  cap  and  the  cryostat  specimen  post  provided 

thermal  contact  between  these  two  components.  Thermal 

contact  between  the  top  cap  and  the  specimen  is  established 

by  the  helium  gas  inside  the  "Kovar"  to  glass  seal.  We 

have  found  that  thermal  contact  is  somewhat  improved  at 

4 

temperatures  below  the  lambda  point  of  He  .  The  reason 

for  this  is  that  at  temperatures  above  the  lambda  point 

thermal  contact  is  maintained  by  gas  conduction,  while  at 

4 

lower  temperatures  a  film  of  superfluid  He  acts  as  the 


thermal  link. 
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FIGURE  5 . 6 


Schematic  showing  the  method  of  specimen  mounting. 
This  drawing  is  approximately  two  times  actual  size. 
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C.  Procedure 

The  operation  and  cool-down  procedure  of  the  cryo¬ 
stat  has  already  been  described  in  chapter  IV.  The  speci¬ 
men  resistance  as  well  as  one  resistance  thermometer 
calibration  point  was  obtained  at  4°K  before  the  exchange 
gas  was  pumped  out  of  the  experimental  chamber  of  the 

cryostat.  After  pumping  out  the  exchange  gas  for  about 

-6 

one  hour  (to  a  gauge  pressure  of  about  10  mm  of  Hg.),  the 
apparatus  was  cooled  to  a  temperature  below  the  super¬ 
conducting  transition  temperature  of  the  specimen.  This 

a 

was  done  with  most  the  He  in  the  cold  part  of  the  cryostat. 

7  7  4 

The  He  gas  in  the  pumping  tube  between  the  He  and  He 

chambers  acted  as  a  thermal  link.  Somewhere  between  4°K 
and  2°K  a  second  resistance  thermometer  calibration  point 
was  obtained.  The  third  resistance  thermometer  calibra¬ 
tion  point  was  normally  obtained  at  about  1.25°K.  If 
temperatures  below  1.25°K  were  required,  the  temperature 
was  further  reduced  to  0.5°K  by  pumping  on  the  He-V 
Measurements  on  a  specimen  were  first  carried  out  with  it 
in  the  superconducting  state.  Measurements  with  the  speci¬ 
men  in  the  normal  state  were  obtained  at  the  end  of  the  run, 
at  4°K,  with  exchange  gas  in  the  cryostat.  Since  normal 
state  measurements  were  only  required  for  relative  conduc¬ 
tance  measurements  this  procedure  enabled  us  to  examine 
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the  oscilloscope  photographs  of  6V-V  curves  obtained  with 

the  specimen  in  the  superconducting  state  and  then  to 

take  photographs  of  oscilloscope  traces  with  the  specimen 

in  the  normal  state  using  the  same  scale.  Thus  the  scale 

used  was  determined  by  the  finest  structure  visible  on 

6V-V  in  the  superconducting  state.  At  the  end  of  the  run., 

with  exchange  gas  in  the  experimental  chamber  of  the  cryo- 

stat  all  the  He^  was  pumped  back  into  storage  and  all  the 

valves  of  the  He  system  were  closed.  The  warm-up  procedure 

4 

from  then  on  was  the  same  as  that  for  any  He  cryostat. 
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CHAPTER  VI 


EXPERIMENTAL  RESULTS 

A.  Current-Voltage  Curves 
(1)  Pb-Al203-Al: 

This  type  of  junction  was  the  first  to  be  attempted. 
Since  the  transition  temperature  of  lead  is  7.2°K  we  were  able 
to  familiarize  ourselves  with  the  method  and  with  the  char¬ 
acteristics  of  good  specimens  by  immersing  the  junction  in 
liquid  helium  in  a  storage  dewar .  For  these  immersion 
experiments  the  specimen  was  mounted  in  a  stainless  steel 
tube  filled  with  He^  gas  at  low  pressure.  Specimens 
cooled  without  such  mounting  were  invariably  ruined  even 
when  dipped  into  liquid  air.  The  sudden  thermal  shock 
that  occurs  when  the  unmounted  specimen  is  dipped  creates 
a  rapid  differential  contraction  between  the  metal  film  and 
the  glass  substrate  which  breaks  the  film.  When  we  developed 
our  experimental  technique  sufficiently,  a  junction  was 

mounted  as  described  in  the  previous  chapter  and  soldered 

4 

to  the  outside  surface  of  the  He  chamber.  The  I-V  curves 
obtained  with  this  specimen  are  shown  in  Fig.  6.1.  The 
specimen  resistance,  recorded  in  the  normal  state,  during 
warm-up  was  about  4000  ohms.  The  room  temperature  resis¬ 
tance  obtained  preceeding  the  run  was  about  1250  ohms. 


. 
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FIGURE  6.1 


Current-Voltage  curves  for  a  Pt>- Al^O^- A1 
junction . 
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This  change  in  resistance  arises  because  more  than  one 
conduction  mechanism  across  the  barrier  is  present. 

One  conduction  mechanism  is  provided  by  electrons  either 
in  the  valence  band  or  impurity  levels  of  Al^O^  that 
are  thermally  excited  into  the  conduction  band  at  room 
temperature  thus  enhancing  the  conductance  of  the  tunnel 
junction.  This  mode  of  charge  transfer  is  greatly  reduced 
at  liquid  air  or  lower  temperatures  since  the  thermal 
energy  available  for  excitation  is  much  lower.  The  second 
conduction  mechanism,  provided  by  tunneling  through  the 
barrier,  is  almost  temperature  independent.  A  thorough 
discussion  of  electron  tunneling  through  thin  insulating 
films  is  given  by  Holm  (1951)  and  by  Fisher  and  G-iaever 
(1961),  In  most  of  our  specimens  the  resistance  at  liquid 
air  temperatures  was  two  to  four  times  the  room  temperature 
resistance,  while  there  was  a  further  increase  in  resistance 
of  only  about  ten  percent  between  liquid  air  and  liquid 
helium  temperatures. 

An  estimate  of  the  thickness  of  the  insulating 
barrier  was  also  made.  It  was  assumed  that  the  equivalent 
circuit  of  the  tunnel  junction  could  be  represented  by  a 
resistor  in  parallel  with  a  parallel  plate  capacitor,  the 
distance  between  t he  plates  being  the  thickness  of  the 
oxide  layer.  The  data  for  this  calculation  were  obtained 
from  the  I-V  curves  at  a  number  of  frequencies.  These 
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measurements  showed  a  capacity  gf  about  0.24pF.  Using  the 
value  of  6.4  obtained  by  Feldman  and  Hacskaylo  (1962)  for 
the  dielectric  constant  of  Al^O^  we  obtain  an  estimate  of 

O 

about  20A  for  the  barrier  thickness.  It  must  be  pointed 
out  that  this  estimate  is  probably  extremely  crude,,  because 
of  the  following  major  uncertainties : 

(i)  Measurements  of  dielectric  constants  of 
films  of  thicknesses  comparable  to  these  would  be  difficult 
to  make;  in  fact  they  would  involve  a  parallel  plate 
capacitor  arrangement ,  not  unlike  our  tunnel  junctions  and 
the  determination  of  film  thickness  would  be  uncertain. 

The  value  of  dielectric  constant  used  has  been  obtained 

by  Feldman  and  Hacskaylo  on  material  which  was  100|~i  thick, 
they  also  pointed  out  that  extrapolation  to  thinner  films 
is  uncertain. 

(ii)  Since  capacitance  measurements  yield  some 
average  thickness  for  the  barrier,  while  tunneling  is 
exponentially  dependent  on  the  barrier  thickness,  any 
width  infered  from  capacity  measurements  is  likely  to  be 
larger  than  the  effective  tunneling  width.  It  is  probably 
safe  to  conclude  that  the  thickness  of  the  AlpO^  layer 

is  less  than  30  A.  This  takes  into  account  the  possibility 
of  having  underestimated  the  dielectric  constant  by  50$. 

Finally  it  should  be  mentioned  that  the  resistance 
of  the  tunnel  junctions  was  found  to  increase  with  time  at 
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room  temperature „  This  increase  was  fairly  rapid  during 
the  first  few  hours  after  the  junction  was  made.  For 
example.,  one  newly  made  junction  with  a  resistance  of 
10  ohms  had  a  resistance  between  100  and  1000  ohms  a 
few  hours  later.  The  rate  of  increase  of  resistance 
decreased  with  time.,  but  no  actual  saturation  was  observed. 
This  type  of  behaviour  occured  in  all  our  specimens  to  a 
greater  or  lesser  degree.,  not  just  in  Pb-AlpO^- A1  junctions. 
If  the  specimens  were  kept  at  liquid  air  or  lower  temper¬ 
atures  their  resistance  remained  constant.  Storing  the 
specimens  in  an  inert  atmosphere  (helium  gas)  at  room 
temperature  did  not  inhibit  the  resistance  increase.  It 
seems  evident  that  whatever  process  caused  the  resistance 
to  increase  depends  on  the  amount  of  thermal  energy  avail¬ 
able.  Assuming  that  the  increase  in.  resistance  is  due  to 
a  thickening  of  the  barrier  layer,,  two  tentative  explana¬ 
tions  are  offered: 

(i)  Oxygen  trapped  between  the  two  metal  films 
combined  with  the  metal  substrate  to  form  more  oxide, 
thus  increasing  the  barrier  thickness. 

(ii)  Assuming  that  the  interface  between  the 
metals  and  the  oxide  layer  is  initially  not  very  sharp, 

a  diffusion  process  in  which  the  metal  moves  in  one  direc¬ 
tion  while  the  oxide  moves  in  the  other  in  such  a  way  as 
to  sha.rpen  the  boundary  would  also  have  the  effect  of 
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increasing  resistance. 

One  or  both  of  the  above  explanations  are  possible 
because  each  of  these  mechanisms  would,  slow  down  as  the 
temperature  is  lowered. 

(2)  Zn-Al203-Al: 

Measurements  on  this  type  of  junction  were  carried 
out  since  no  measurements  of  the  energy  gap  in  the  density 
of  electron  states  could  be  found  in  the  literature  for 
superconducting  zinc.  The  I-Y  characteristics  obtained  at 
0.5°K  are  shown  in  Fig,  6.2. 

(3)  Similar  measurements  were  also  carried  out  for  aluminum 
and  indium  using  Al-A^O^-Al,  (see  Fig.  6.3)  and  In-Al^Cy-Al 
tunnel  junctions . 

B.  The  Energy  Gap 

Values  of  the  energy  gap  in  the  density  of  elec¬ 
tronic  states  of  a  superconductor  can  be  obtained  from  the 
current-voltage  curves.  For  combinations  of  a  normal 
metal  and  a  superconductor  the  voltage,  =  A/e,  corres¬ 
ponds  to  the  point  at  which  the  relative  conductance  is 
unity.  In  the  case  of  two  dissimilar  superconductors  the 
voltage,  V2  =  (A1  +  A0)/e,  is  given  by  the  point  where  the 
I-V  curve  is  steepest,  while  the  sharp  peak  preceeding  the 
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FIGURE 
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Current-Voltage  curves  for  a  Zn-Al^O^-Ai  tunnel 
junction  at  T  =  0„5°K. 
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FIGURE  6.3 


Current- voltage  curves  for  an  Al-Al^O.-,-  A1  tunnel 
junction  at  T  =  0o5°K. 
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negative  resistance  region  occurs  at  a  voltage,, 

=  (Ax  -  Ag)/e.  Values  of  the  energy  gap  so  obtained 
are  tabulated  below. 


TABLE  I 


Values 

of  the 

energy 

gap  (2A)  in  units 

T  =  0.5°K 

of  (10  ^eV)  at 

Metal 

This 

Work 

Calculated  from 

Giaever 

and 

the  BCS  theory* 

Megerle 

(1961) 

Pb** 

2.47 

-  0.10 

to 

IV) 

0 

2.68  + 

0 

0 

ov 

In 

1 — 1 

b 

VO 

-0,05 

1.03 

1.05  - 

0.03 

A1 

0.39 

±  0.03 

0.355 

0.42  - 

0.06 

Zn 

0.23 

-  0  0  02 

0.234 

—  —  — 

-  —  - 

C.  Measurements  of  Relative  Conductance 

Giaever,  Hart.,  and  Megerle  (1962)  and  Rowell, 
Chynoweth,  and  Phillips  (1962)  have  reported  structure 
in  the  relative  conductance  curve  of  lead.  The  observed 


*  In  these  calculations  the  transition  temperatures 
given  by  Cochran  and  Mapother  (1961)  were  used. 

**  The  gap  energy  in  lead  was  obtained,  at  1.24°K  but  is 
not  expected  to  change  significantly  at  lower  temper¬ 
atures  because  of  the  high  transition  temperature  of 
lead . 
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structure  was  found,  to  be  harmonic.,  occuring  at  energies 
A  +  n  E^.,  where  E^_  is  the  energy  of  the  transverse  acoustic 
phonon  in  lead  and  n  is  an  integer .  Since  this  structure 
has  until  now  only  been  observed  in  lead  it  is  of  interest 
to  examine  other  superconductors.  Since  such  multiphonon 
effects  should  be  easiest  to  observe  in  superconductors 
having  a  strong  electron-phonon  interaction  the  most  logical 
metal  to  investigate  is  mercury.  Since  we  were  not  success¬ 
ful  in  making  tunnel  junctions  with  mercury  we  decided  to 
investigate  indium  and  aluminum,  in  which  the  electron- 
phonon  interaction  is  not  quite  as  strong. 

(l)  Indiums* 

The  relative  conductance  of  indium  was  measured 
using  an  In-Al^O^-Al  tunnel  junction.  Specimens  used  had 
tunneling  resistances  of  30,  l4o,  and  230  ohms.  A  small 
depression  appears  (between  -12  and  -l6  mV)  inside  the 
dotted  square  in  Fig.  6.4  (a).  This  corresponds  to  struc¬ 
ture  in  the  relative  conductance  of  indium.  There  is  an 
exactly  similar  effect  occuring  in  the  positive  half  of 
the  curve,  but  is  not  as  readily  visible  due  to  the  steep¬ 
ness  of  the  dl-V  curve  there.  The  asymmetry  in  the  6V-V 
curves  occurs  in  all  tunnel  junctions  of  the  type  A-Ax0y-B, 


See  also  Appendix  IV 
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FIGURE  6.4 

Relative  conductance  of  an  In-AlgO^-Al  tunnel  junction. 

(a)  6Y-V  traced  from  an  oscilloscope  photograph. 

This  photograph  was  taken  with  small  sensitivity 
in  order  to  show  a  large  voltage  range. 

(b)  Relative  conductance  of  a  230  ohm  specimen  at 

T  =  2.04°K.  Indium  but  not  the  aluminum  is  super¬ 
conducting  . 

(c)  Relative  conductance  of  a  30  ohm  specimen.  Both 
indium  and  aluminum  are  superconducting. 
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where  A  and  B  are  dissimilar  metals.,  and  probably  arises 
because  of  the  difference  in  work  function  of  the  two 
metals.  The  same  asymmetry  appears  in  both  the  normal 
and  superconducting  states  and  hence  does  not  affect  the 
results.  The  relative  conductance  curves  Pig.  6.4  (b) 
and  (c)  were  evaluated  from  6V-V  curves  using  the  sensi¬ 
tivity  of  1  pV/cm  for  6V.  The  relative  conductance  is 
then  determined  to  an  accuracy  of  -  0.0001.  The  level  of 
the  whole  relative  conductance  curve  with  respect  to 
unity  can  be  ascertained  only  within  4  0.0005^  because 
the  steepness  of  the  6V-V  curves  at  the  higher  voltages 

makes  accurate  measurement  of  6V  -  6V  difficult.  The 

s  n 

dashed  line  is  the  theoretical  curve  calculated  from  the 
BCS  theory  for  T  =  0°K.  The  curve  seems  to  remain  asymptotic 
to  unity  for  voltages  greater  than  =  20  mV,  but  above 
about  30  mV  heating  effects  become  noticeable.  Heating 
effects  cause  discrepancies  between  the  trace  obtained 
with  increasing  sweep  voltage  and  that  obtained  with 
decreasing  sweep  voltage.  If  heating  effects  are  present 
then  the  slower  the  sweep  the  larger  the  discrepancy.  Of 
the  specimens  used  none  showed  heating  effects  below  =  20  mV. 
The  specimen  with  the  highest  tunneling  resistance  (Pig.  6.4 
(b))  was  the  least  affected  by  heating,  this  is  because  the 
power  input  to  the  junction  varies  inversely  with  junction 


resistance . 
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FIGURE  6.5 


Relative  conductance  of  an  In-AlpO^-Al  tunnel  junction 
in  the  vicinity  of  the  energy  gap. 
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We  also  examined  the  relative  conductance  in  the 
vicinity  of  the  energy  gap.  Our  results  are  in  agreement 
with  those  of  Giaever,  Hart,  and  Megerle  (1962)  and  are 
shown  in  Pig.  6.5.  Because  of  the  steepness  of  the 
relative  conductance  curve  in  this  region  we  cannot  obtain 
as  great  an  accuracy  at  these  small  biases  as  we  can  at 
higher  voltages  (^6  to  ^20  mV) .  Nevertheless  the  data  in 
Fig.  6.5  are  accurate  to  about  -  5$  for  voltages  below 
^  1  mV,  and  becomes  more  accurate  as  the  voltage  across 
the  tunnel  junction  increases. 

( 2 )  Aluminum : 

Relative  conductance  measurements  were  made  on 
Al-Al^O^-Al  tunnel  junctions,  a  portion  of  one  of  the 
resulting  curves  is  shown  in  Pig.  6.6.  Since  we  wanted 
to  carry  our  measurements  out  to  voltages  which  correspond 
to  energies  higher  than  the  Debye  energy  of  aluminum, 
tunnel  junctions  were  made  having  high  resistances  (200 
to  2000  ohms)  to  minimize  heating  effects*.  As  can  be 
seen  from  Fig.  6.6  agreement  with  the  BCS  theory  is 
surprisingly  good.  (The  dashed  line  is  the  relative  con¬ 
ductance  calculated  using  equation  (A-6)).  Small  humps 
and  depressions  in  the  relative  conductance  appear 
throughout  the  whole  curve  out  to  about  60  mV.  The 
amplitude  of  these  effects  varies  starting  from  a  rather 

*  The  Debye  temperature  of  aluminum  is  430°K  which 
corresponds  to  a  voltage  of  37  mV. 
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FIGURE  6.6 

an  Al-Al^O^-Al  tunnel 

curve  indicates  the 
theory  at  the  absolute 


Relative  conductance  curve  for 

junction  at  T  =  0.5°K.  Dashed 
behaviour  predicted  by  the  BCS 
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small  amplitude  near  the  gap  region  and  increasing  some¬ 
what  at  larger  biases.  Several  attempts  to  find  significant 
regularities  or  periodicities  in  the  energies  at  which  the 
humps  occurred  were  unsuccessful.  These  humps  were  clearly 
visible  on  the  6V-V  curves  and  occured  only  when  the 
aluminum  was  in  the  superconducting  state.  They  were  quite 
reproducible  and  several  photographs  showed  exactly  the 
same  structure  on  separate  occasions,  as  long  as  the 
specimen  was  kept  at  liquid  helium  temperatures.  The 
same  humps  persisted  on  different  specimens,  but  humps 
having  the  same  character  appeared  at  higher  energies 
(never  more  than  10  3  eV)  in  a  specimen  that  was  allowed 
to  warm  up  to  room  temperatures  between  runs. 


■  :  .  ■ 
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CHAPTER  VII 

CONCLUSIONS  AND  DISCUSSION 

A.  The  Energy  Gaps 

Table  I  shows  that  the  values  obtained  for  the 
energy  gap  in  the  density  of  electron  states  are  in 
remarkable  accord  with  the  values  predicted  by  the  BCS 
theory.  This  good  agreement  with  the  BCS  theory  seems 
characteristic  of  thin  film  experiments,  such  as  those  of 
Giaever  and  Megerle  (1961),,  Giaever,  Hart  and  Megerle 
(1962)  *  Townsend  and  Sutton,,  (1962).,  and  of  Ginsberg  and 
Tinkham  (i960).  On  the  other  hand  measurements  of  the 
energy  gap  in  single  crystals,,  done  by  e.g.  Morse,  Olsen, 
and  Gavenda  (1959) *  Morse  and  Bohm  (i960).  Mackintosh 
(i960),  using  ultrasonic  attenuation  indicate  dependence 
of  the  energy  gap  on  crystalline  orientation.  Such 
anisotropy  is  not  surprising  when  one  considers  the  com¬ 
plex  Fermi  surfaces  of  metals  which  become  superconductors. 
It  may  be  that  anisotropy  does  not  manifest  itself  very 
strongly  in  thin  film  experiments  because  of  film  thick¬ 
nesses  that  are  less  than  the  coherence  length.  In  such  a 
case  the  wave  functions  for  the  electrons  in  the  super¬ 
conducting  state  might  be  affected  to  an  extent  large 
enough  to  mask  any  anisotropy  effects.  Such  a  geometrical 
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restriction  of  the  electronic  wave  functions  might  be 
expected,  to  modify  some  of  the  parameters  which  are 
characteristic  of  the  superconducting  state.,  such  as  the 
energy  gap . 

As  can  be  seen  from  Table  I,  our  results  for 
the  energy  gap  of  lead  are  about  12$  higher  than  those 
predicted  by  the  BCS  theory.  This  is  because  the  theory 
is  only  applicable  in  the  weak  coupling  limit.,  that  is 
for  metals  in  which  the  electron-phonon  interaction  is 
moderately  weak.  A  measure  of  the  strength  of  this  inter¬ 
action  is  usually  expressed  in  terms  of  the  coupling 
constant,  Tc/^]y  where  0^  is  the  Debye  temperature.  The 
weak  coupling  limit  of  the  BCS  theory  implies  T  £  0.01, 
whereas  for  lead  Tc/0^  =  0.075. 

For  each  of  our  aluminum  tunnel  junctions  we 
found  that  the  superconducting  transition  temperature  was 
higher  than  that  of  bulk  aluminum.  This  is  in  accord 
with  the  observations  of  G-iaever,  Hart,  and  Megerle  (1962) 
which  showed  that  the  transition  temperature  increased  as 
the  film  thickness  decreased.  These  workers  also  noted  that 
the  energy  gaps  obtained  from  their  tunneling  measurements 
were  somewhat  higher  than  those  predicted  by  the  BCS  theory 

-X* 

for  aluminum. '  They  attributed  this  behaviour  to  the 
effect  of  internal  stresses  on  their  films.  Our  results 
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For  aluminum  Tc/0p  =  0.0027.,  thus  the  weak  coupling 
limit  of  the  BCS  theory  should  apply. 
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(Table  I)  are  In  general  accord  with  those  of  Giaever 
and  his  co-workers^  but  we  feel  that  a  more  satisfactory 
explanation  of  the  high  transition  temperatures  might  be 
found  in  terms  of  the  coherence  length.  In  the  BCS 
theory  the  relation  for  the  coherence  length,  ^  in 
terms  of  the  energy  gap  parameter,  A  ,  is 


hv. 


irAr 


(7.1) 


where  v  is  the  Fermi  velocity.  Since  these  films  all 
have  a  thickness  which  is  less  than  the  coherence  length, 
the  bulk  coherence  length  might  be  replaced  by  a  new 
coherence  length  £,(t)  <  where  t  is  the  film  thickness. 

The  transition  temperature  would  then  still  be  linked  to 
the  energy  gap  by  the  BCS  relation  2A0  =  3-52  kTc .  Both 
indium  and  zinc  show  extremely  good  agreement  with  the 
energy  gap  value  predicted  by  the  BCS  theory,  indeed  so 
good  that  in  the  light  of  what  has  been  said  above  it 
may  be  somewhat  fortuitous. 


B.  The  Relative  Conductance  Measurements 

A  discussion  of  the  relative  conductance  of 
indium,  which  has  been  accepted  for  publication,  is 
given  in  Appendix  IV,  and  only  the  results  for  aluminum 
will  be  discussed  here.  A  portion  of  the  relative  con- 
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ductance  curve  for  an  Al-Al^O^-Al  tunnel  junction  Is 
shown  in  Fig.  6.6.  The  agreement  between  the  theoretical 
curve  and  the  experimental  points  is  extremely  close., 
being  better  than  one  percent  throughout.  Measurements 
were  made  to  biases  of  about  60  mV,  but  only  a  portion  of 
the  results  were  plotted  as  the  remainder  are  asymptotic 
to  unity  and  show  about  the  same  scatter.  As  mentioned 
above  (Chapter  VII,  Section  C)  we  were  unable  to  find  any 
regularity  in  the  observed  structure.  It  is  at  first 
sight  tempting  to  attribute  this  structure  to  multiphonon 
effects  similar  to  those  observed  by  Rowell,  Chynoweth,  and 
Phillips  (1962)*,  but  which  appear  at  various  energies 
due  to  the  real  lattice  spectrum  not  conforming  to  an 
Einstein  approximation.  However  if  these  were  in  fact 
multiphonon  effects  one  would  expect  the  amplitude  of  the 
peaks  to  diminish  as  the  electron  energy  increased,  particu¬ 
larly  above  the  Debye  energy,  since  the  probability  of  a 
two  phonon  process  is  less  than  that  of  a  single  phonon 
process.  An  explanation  for  these  small  bumps  in  the 
relative  conductance  curve  is  not  apparent. 

C.  Suggestions  for  Further  Investigation 

The  above  investigations  in  this  relatively 
new  field  of  electron  tunneling  into  superconductors  are 
naturally  of  somewhat  preliminary  nature.  The  experience 


* 


See  Chapter  II,  Section  D. 
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gained  thus  far  has  however  clarified  the  scope  of  the 
technique  and  stimulated  ideas  for  further  experiments  that 
should  provide  definitive  information.  Some  of  these  are 
hriefly  listed  below: 

(i)  The  measurement  of  energy  gaps  by  electron 
tunneling  into  bulk  superconductors ,  using  single  crystals 
of  known  orientation.  Such  measurements  should  reveal  the 
anisotropy  of  the  energy  gap  and  its  relation  to  the  Fermi 
surface  of  the  metal  in  its  normal  state. 

(ii)  Electron  tunneling  experiments  with  the 
transition  metal  superconductors  ruthenium  and  osmium. 

Such  experiments  are  of  interest  from  a  theoretical  stand¬ 
point  because  these  metals  show  no  isotope  effect,  which 
is  an  indication  that  superconductivity  in  them  does 

not  arise  from  the  electron-phonon  interaction. 

(iii)  Electron  tunneling  experiments  with  super¬ 
conducting  alloys.  These  experiments  would  provide  infor¬ 
mation  about  the  effects  of  alloying  on  the  attractive 
interaction  which  gives  rise  to  superconductivity. 

(iv)  Measurements  of  the  energy  gap  in  thin 
films  along  with  a  careful  determination  of  both  the 
transition  temperature  and  film  thickness,  to  establish 
the  correlation  between  these  parameters. 

(v)  A  careful  re-examination  of  the  relative 
conductance  of  lead  to  determine  whether  any  further 
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multiphonon  structure  besides  that  observed  by  Rowell  et 
al .  exists.  Neither  this  group,,  nor  Giaever  and  his  co¬ 
workers  have  indicated  the  sensitivity  with  which  they 
are  able  to  observe  the  relative  conductance,  and  it  is 
quite  possible  that  our  bridge  technique  is  an  order  of 
magnitude  more  sensitive.  Furthermore  efforts  should  be 
made  to  produce  tunnel  junctions  in  which  one  component 
is  mercury  and  determine  if  it  exhibits  multiphonon  effects 
as  well. 

D.  A  Proposed  New  Application  of  Electron  Tunneling  to 
Normal  Metals 

The  value  of  electron  tunneling  techniques  as 
a  solid  state  research  tool  stems  from  the  fact  that  by 
virtue  of  the  barrier  layer  one  can  inject  electrons  of  a 
known  energy  and  direction  of  motion  into  a  metal.  Con¬ 
sider  a  tunnel  junction  of  the  type  A/  barrier/  A',  where 
A  is  a  single  crystal  of  known  orientation  and  A'  is  an 
evaporated  film  of  metal  A.  If  we  apply  a  positive  volt¬ 
age  to  metal  A,  then  only  electrons  with  a  direction  of 
travel  that  is  essentially  perpendicular  to  the  barrier 
layer  will  tunnel  from  A1  to  A  because  of  the  exponential 
dependence  of  the  tunneling  probability  on  the  barrier 
thickness.  Thus  the  direction  of  motion  with  respect  to 
the  crystalline  axes  of  metal  A  is  known.  The  wave-vector 
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of  the  electron  is  now  fully  determined,  and  if  it  corre¬ 
sponds  in  energy  to  the  Brillouin  zone  boundary  for  this 
wave  vector,  the  electron  will  suffer  a  Bragg  reflection, 
that  is  it  will  not  be  permitted  to  tunnel  at  all .  In  a 
practical  situation  there  will  be  some  contribution  to  the 
current  even  at  these  energies  for  two  reasons.  (i)  The 
electron  may  be  scattered  by  impurities  or  lattice  imper¬ 
fections.  (ii)  The  electron  may  be  scattered  by  lattice 
vibrations.  For  this  reason  such  experiments  must  be 
carried  out  at  low  temperatures  (liquid  helium)  where  the 
phonon  density  is  low.  Also  the  single  crystal  of  mater¬ 
ial  A  must  be  very  pure  and  have  a  very  carefully  prepared 
surface  so  that  lattice  irregularities  scatter  only  a 
small  fraction  of  the  tunneling  electrons.  Such  an  experi¬ 
ment  enables  one  to  determine  the  energies  involved  in 
umklapp  processes  (Bragg  reflections)  directly.  This 
information  would  be  of  use  in  determining  the  umklapp 
contribution  when  calcuating  transport  coefficients.  If 
such  measurements  are  carried  out  in  many  crystal  direc¬ 
tions  and  if  one  assumes  a  general  form  for  the  dispersion 
relation  (i.e.  energy  as  a  function  of  wave  vector)  one 
can  in  principle  determine  the  Fermi  surface. 
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APPENDIX  I 


EVALUATION  OF  SOME  INTEGRALS  PC CUR I NO  IN 
TUNNELING  PROBLEMS 


Evaluation  of  I  =  aj  Q( E) pn ( E) p0 (E+eV) [ f (E) -f ( E+eV) ] dE, 
at  T  =  0°K  for  small  applied  voltage,,  V.  and  arbitrary  Q(E)  . 


At  T  =  0°K 


f(E)  =  0 
f(E)  =  1 
f ( E+e  V )  =  0 
f ( E+eV)  =  1 


for  E  >  0 
for  E  <  0 
for  E  >  -eV 
for  E  <  -eV 


Without  loss  of  generality  V  is  assumed  positive  throughout 

o  -eV 

.  .  I  =  af  Q(E)Pl(E)pQ(E+eV)dE  -  a  Q( E)  pn  ( E)  Pp  ( E+eV)  dE 


-  00 
O 


i 


a 


/_ev^( E) Pi  (E)pQ(E+eV)dE. 


Carrying  out  a  Taylor  expansion  about  the  Fermi  level: 


I  =  a/7vQ(E)[ppo)+E(3Pl/ciE)E=0+.  ..][p2(eV)+E(3p2/3E)E=eV+.  .  .  ]  dE 

=  a/°  Q(E)p,  (0)p„(eV)dE  +  .... 

-e  V  1  ^ 

Assuming  the  densities  of  states  (p^  and  p^)  to  be  approxi¬ 
mately  constant  in  the  vicinity  of  the  Fermi  surface,  higher 
order  terms  may  be  neglected,  furthermore  if  eV«£p,  where 
is  the  Fermi  energy,  the  integral  reduces  to: 

I1 


I K  I 


108 . 


o 

I  -  ap1(0)p2(0)  /eV Q(e)  dE.  (A-l) 

(i)  For  both  metals  In  the  normal  state  we  set  Q(E)  =  1 
hence  Inn  -  ap1(0)p2(0)  eV.  (A-2) 

(li)  For  the  case  In  which  one  of  the  metals  Is  In  the 
normal  state  while  the  other  Is  in  the  superconducting 
state,  one  obtains  on  the  basis  of  the  BCS*  theory: 

Q(E)  =  RP  |  E  |  (E2~A2)-1/2 

J°eV  Q(E)dE  =  RP  j£v|E  |  (E2-A2)"l/2dE 
=  RP  /^E(E2-A2)_1/2dE 
=  RP[  (eV)2-A2]  1/7 ^ 

Substituting  into  equation  (A-l)  we  have 

Ins  =  0  for  eV  <  A  (A-3) 

Ins  =  ap1(0)p2(0)[(eV)2-A2]1/2  for  eV  >  A  (A- 3a) 


(iii)  For  the  case  of  two  identical  superconductors,  at 
T  =  0°K,  we  again  assume  the  BCS  theory. 


Q(E) 


p^p  _ I  E  |  |  E  +  eV  | _ 

/(E2-A2) [ (E+eV)2  _  A2] 


o 

/  Q(E)dE  =  RP 
-eV 


|-°  _ _ lE[|  E+eV  |  dE 

"eV(E+A)  (E-A)  (E+eV-A)  (E+eV+A) 


=  RP  / 


eV 

o 


_ E  (eV-E) _ 

■/(-E+A) ( -E-A) (-E+eV-A) (-E+eV+A) 


*  This  is  the  original  BCS  theory,  Bardeen,  Cooper,  and 
Schrieffer  (1957) ^  in  which  the  energy  gap  parameter.  A, 
is  independent  of  energy. 
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o 

/  Tr  Q(E)dE  =  0  for  eV  <  2A  because  the  denominator  becomes 
J  -eV 

imaginary.  For  eV  >  2A  we  have: 
e  V  pV 

J  =  /  Q(E)dE  =  RPJ  - _ -  ■Ul/g.ey-E),  1/geV]  dE  ,  „  ,  _ _ _ 

0  °  Vi  (l/2eV-E)  +  (A-l/2eV)  ]  [  ( l/2eV-E)  -  ( A+l/2eV)  ] 

.  j  e  V  _ [l/2eV-(l/2eV-E)  ]dE  _ _ 

°  Vt  l/2eV-E)  +  (A+l/2eV)  ]  [  ( l/2eV-E)  -  ( A-l/2eV)  ] 


put:  u  =  l/2eV  -  E 

a  =  1/2 (2A  -  eV) 
p  =  1/2 (2A  +  eV) 


then : 


J 


RP  / 


1/2eV  [(l/2eV)g-ug1  du 

V2eV  i/(u+a)  (u-p)  (u+p) (u-a) 


RP  f 


l/2eV 

l/2eV 


[ (l/2eV)2~u2]du 

y(u2-a2)  (u2-p2) 


RP-1-  f1//2eV  [  (l/2eV)2-u2]du 

af3  -l/2eV  V (u2/a2-  l) (u2/p2  _  X) 

a2  1/2fcV  [ ( eV/2a) 2  -  (u/a) 2] d(u/a) 
p2  -l/2eV  -v/(  i-u2/d2)  ( l-a2/p2  •  u2/a2) 


put:  t  =  u/a  and  -k  =  a/p  =  (2A  -  eV)/(2A  +  eV) 


J  = 


/  p  a  2  eV-'^A 
RP  1/2  i2AAg.Al-  f 


p 

t  at 


(2A+eV)  V'(l-t2)(l-k2t2' 


+ 


(eV) 


eV 


dt 


•eV  / 


v/(l-t2) ( l-k2t2) 


eV+2A 


eV-2A 
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Since  eV  >  2A  and  both  A  and  eV  are  real  and  positive,,  the 
real  part  of  the  integral  has  contributions  only  in  the 
range  0  ^  t  <  1,  hence: 


J  -  (2A+eV)k2  f1 


t2dt 


+  (eV)2  f 1 


dt 


o  v/(l-t2) (i_k2t2)  eV+2A  °  v/( 1-t2) ( l-k^t2) 

using  the  relations:*  l/k2[K(k)  -  E(k)]  =  /' 


■  1  t^dt 


v/(l-t2)(i-k2t2) 


and  K(k)  =  f 


1 


dt 


°v/(l-t2)  (i_k2t2) 
were  K  and  E  are  complete  elliptic  integrals  of  the  first 
and  second  kind  respectively. 

Thus  : 


j  =  ,_(_g.Y)_.  k(c-_v:  j-)  +  (eV+2A)  [E(-T-f— )  -  K(^:^)  ] 


eV-2A' 


eV-2A' 


eV-2A' 


eV+2A  eV+2A 


eV+2A' 


eV+2A' 


I  0  =  0  for  eV  <  2A 
ss 


(A-4) 


Iss  aP 


(0) 


leXlf  K(eV-2A} 
eV+2A  veV+2A; 


+  (eV+2A)[E(|^|$)-K(|^||)]  eV>2A 


eV+2A' 


( A-4a) 


which  agree  with  those  obtained  by  Giaever  Hart,,  and 
Megerle  (1962). 


The  evaluation  of  the  I-V  characteristics  for 
non- zero  temperatures  are  much  more  difficult  and  can 
only  be  solved  by  numerical  methods.  Shapiro.,  Smith,  Nicol, 
Miles,  and  Strong  (1962)  discuss  such  numerical  integrations. 


* 


Magnus  and  Oberhettinger  (1954),  pp .  112-3 


’i£ 
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APPENDIX  II 

v 

NUMERICAL  EVALUATION  OF  THE  NORMALIZED  CURRENT 
AND  RELATIVE  CONDUCTANCE  FOR  THE  CASE  OF  TWO 
IDENTICAL  SUPERCONDUCTORS  AT  THE  ABSOLUTE  ZERO. 


Defining  the  normalized  current,,  J  =  I  /I  , 

S  S  S  o'  liri 

and  the  relative  conductance  Cgs  =  (dl/dV) ss/(dl/dV)nn„  for 
the  case  of  two  identical  superconductors  equation  (A-4a) 
becomes : 

J  =  K(k)  +  E(k)  _  K(k)  (A-5) 

ss  x+2  x  v  ' 


G 


ss 


x(x+4) 

(x+2)2 

+  (x+2) 


K(k)  +  [Uk  .  (x+2)  ] 

x+2  dx 

dEW.  +  E(k)  -  K(k)  . 
dx 


(A- 6) 


Where  x  =  eV/A  and  k  =  (x-2)/(x+2) . 

K(k)  and  E(k)  may  be  evaluated  numerically  using  the  approxi¬ 
mation  formulae  given  by  Hastings  (1955)-  In  order  to 
evaluate  the  relative  conductance  the  numerical  results  of 
(A-5)  may  not  be  differentiated.  The  reason  for  this  being 
that  the  approximation  formulae  used  for  K  and  E  are  in  the 
form  of  Chebyshev  series  which  oscillate  with  small  error 
about  the  true  values.  If  one  puts  y  =  (x-2)  /(x+2) 


then : 


dK 

dy 


1/2 ( y- 1 ) K+E 

y(i-y) 
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dE  _  1/2 (E  -  K) 

dy  y 

may  be  substituted  into  G  ( cf .  Jahnke  and  Emde  1945, 

p.  76)  in  order  to  evaluate  (A-6)  numerically.  The  IBM 

1620  digital  computer  was  used  in  order  to  compile  the 

table  of  JOC!  and  G  given  below, 
s  s  s  s 


TABLE  II 


Energy  in 


I. 


'ss 


( dl/dV) 


ss 


■  s  of  A0  * 

I 

nn 

(dl/dV). 

2.1 

.80376 

1.16424 

2.2 

.81986 

1.15195 

2.3 

.83406 

1.14100 

2.4 

.84664 

1.13119 

2.5 

.85784 

1.12238 

2.6 

.86786 

1.11443 

2.7 

. 87  686 

1.10723 

2.8 

.88497 

1.10069 

2.9 

.89231 

1.09473 

3.0 

.89897 

1.08929 

3.1 

.90503 

1.08430 

3.2 

.91056 

1.07972 

3.3 

.91562 

1.07550 

3.4 

.92026 

1.07161 

3.5 

.92454 

1.06801 

3.6 

.92848 

1.06467 

3.7 

.93212 

1.06158 

3.8 

.935^9 

1.05870 

3.9 

.93862 

1.05602 

A0  is  the 

half  width  of  the 

energy  gap  at  T  = 

nn 


bl 
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Table  II  (Continued) 


4.0 

.94152 

1.05351 

4.1 

.94422 

1.05117 

4.2 

.94674 

1.04898 

4.3 

.94910 

1.04693 

4.4 

•95130 

1.04500 

4.5 

.95336 

1.04319 

4. 6 

.95529 

1.04148 

4.7 

.95711 

1.03988 

4.8 

.95882 

1.03836 

4.9 

.96042 

1.03693 

5.0 

.96194 

1.03558 

5.2 

.96472 

1.03309 

5.4 

.96721 

1.03084 

5.6 

.96944 

1.02882 

5.8 

.97146 

1.02699 

6.0 

.97328 

1.02533 

6.2 

.97493 

1.02381 

6.4 

•97643 

1.02243 

6 . 6 

.97780 

1.02116 

6.8 

.97906 

1.01999 

7.0 

.98021 

1.01892 

7.2 

.98127 

1.01793 

7.4 

.98224 

1.01702 

7.6 

.98314 

1.01617 

7.8 

.98398 

1.01539 

8.0 

.98475 

1.01466 

8.2 

.98547 

1.01398 

8.4 

.98614 

1.01334 

8.6 

.98676 

1.01275 

8.8 

.98734 

1.01220 
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Table  II  (Continued.) 


9.0 

.98789 

1.01168 

9.2 

.98840 

1.01119 

9.4 

. 98888 

1.01073 

9 . 6 

.98933 

1.01030 

9.8 

.98975 

1.00990 

10.0 

.99015 

1.00952 

10.2 

.99052 

1.00916 

10.5 

.99105 

1.00865 

11.0 

.99183 

1.00790 

11.5 

.99251 

1.00724 

12.0 

.99311 

1 . 00666 

12.5 

.99364 

1.00615 

13.0 

.99411 

1.00569 

13.5 

.99453 

1.00529 

14.0 

.99491 

1.00492 

14.5 

.99525 

1.00459 

15.0 

.99556 

1.00430 

15.5 

.99584 

1.00403 

16.0 

.  99609 

1.00378 

16.5 

.99632 

1.00356 

17.0 

.99654 

1.00336 

17.5 

.99673 

1.00317 

18.0 

.99691 

1.00300 

18.5 

.99707 

1.00284 

19.0 

.99722 

1.00270 

19.5 

.99736 

1.00256 

20.0 

•99750 

1.00244 

20.5 

.99762 

1.00232 

21.0 

.99773 

1.00222 

22.0 

.99793 

1.00203 

23.0 

.99811 

1.00186 
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Table  II  (Continued) 


24.0 

.99827 

1.00171 

25.0 

.99841 

1.00158 

26.0 

.99853 

1.00147 

27.0 

.99864 

1.00136 

28.0 

.99874 

1.00127 

29.0 

.99883 

1.00119 

30.0 

.99891 

l.ooill 

31.0 

.99898 

1.00105 

32.0 

.99904 

1.00099 

33.0 

.99910 

1.00093 

34.0 

.99916 

1.00088 

35.0 

.99921 

1.00083 

36.0 

.99925 

1.00079 

37.0 

.99930 

1.00075 

38.0 

.99934 

1.00072 

39.0 

.99937 

1.00068 

40.0 

.99941 

1.00065 

4l.O 

.99944 

1.00062 

42.0 

.99946 

1.00059 

43.0 

.99949 

1.00057 

44.0 

.99952 

1.00055 

45.0 

.99954 

1.00052 

46.0 

.99956 

1.00050 

47.0 

.99958 

1.00048 

48.0 

.99960 

1.00047 

49.0 

.99962 

1.00045 

50.0 

.99963 

1.00043 

51.0 

.99965 

1.00042 

116. 


APPENDIX  III 

MACHINE  DRAWINGS  OF  THE  LOW  TEMPERATURE  PART 
OF  THE  He3  CRYOSTAT 


This  appendix  contains  the  machine  drawings  of 
the  various  parts  of  the  He  cryostat.  These  drawings  are 
those  which  were  submitted  to  the  machine  shop.  It  is  a 
pleasure  to  acknowledge  the  work  of  Mr.  G.  Sanders  in 
machining  all  these  pieces.  Various  ferrules  and  flanges 
used  in  the  room  temperature  part  of  the  apparatus  are  of 
standard  design  and  construction  and  do  not  appear  in  this 
appendix . 


FIGURE  A-l 


Two  sections  showing  assembled  view  of  lower  part 
of  the  He  cryostat* 


SECTION  f  f  SECTION 


k\\\\N 


SS33- 


o 


tf= 

CD 


-scar 


O 


o 


o 


! 


118. 


FIGURE  A- 2 


He 


4 


radiation  trap 


showing  He 


4 


needle  valve  construction. 


This  is  shown  as  A  in  Figure  4.2. 
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FIGURE  A- 3 


3 

HeJ  radiation  trap.  The  two  1/8  inch  ferrules  were 
not  made.  One  of  these  is  to  be  used  in  the  modi¬ 
fication  proposed  in  Chapter  II,  Section  \i,  to 

enable  the  He  condensing  capillary  to  enter  the 
3 

He  pumping  tube.  This  radiation  trap  is  shown  as 
B  in  Fig„  4.2. 
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FIGURE  A- 4 


The  He 


3 


Chamber . 


Shown  as  H  and  J  in  Figure  4.2. 
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FIGURE  A- 5 


Top  of  outer  brass  can,,  shown  as  D  in  Figure  4.2. 
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SECTION  D-  D 


He3  CRYOSTAT 

TOP  OF  OUTER  CAN 
Scale  /> 
tf’A'/A  brass 

J.  Ad  lor  Aug-3,  19(51, 


122. 


FIGURE  A- 6 


The  He  chamber ,  shown  as  F  in  Figure  4.2. 
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FIGURE  A- 7 


Various  ferrules  and  their  location. 


FIGURE  A- 8 


Outer  can.,  shown  as  K  in  Figure  4.2. 
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FIGURE  A- 9 


Two  drawings  showing  the  construction  of  the  dewar 


cap . 


SECTION  I-I 


^ZZZZZZL  d 


A'orfJ: 

(»)  -6cc  4r  _Z.  -  27  t-  o  fc?  c  Ucj,  b»&  f?  i*ess  £:  t°  -5.S .  1> e  27  tG-  . 

(2)  6«t  £■  <  o  *  Hi  -  M  to  be  m  &  c/  c  5  r°  ****  Zy*!'  ^  <■*•+&  3  A  «  «.  C 

it  aro^c^^l  6  Ctf.  fc-ioi^  l'I  a  Mci  Si  l/«r  &  c>  I  of  e  r  C.  c(  ■ 

(3)  *'».*  ioPc  e/ 1  Jt  m  •  t  *  <=>  iec  tie»i  ~!E-  ~  hi  t.o  be  4-  ". 


He3  CRYOSTAT 


DEWAR  CAP 

(drawing  2  ) 

T  ftjlmr-  ^  I. 


126. 


APPENDIX  IV 


SUMMARY  OF  RESULTS  AND  CONCLUSIONS  FOR  THE 
EXPERIMENTS  INVOLVING  INDIUM 


This  appendix  contains  the  manuscript  of  an 
article  submitted  for  publication  in  Physical  Review 
Letters  and  is  now  in  press. 
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Structure  in  the  Density  of  States  of  Superconducting 

Indium  Films*. 

hy 

J.G.  Adler  and  J.S.  Rogers 

Department  of  Physics,  University  of  Alberta 
Edmonton,  Alberta,  Canada 

(Received  l6  January  1963) 


Giaever,  Hart,  and  Megerle1,  and  more  recently 

2 

Rowell,  Chrnoweth,  and  Phillips  have  reported  humps  in 
the  conductance  vs.  voltage  curves  for  electrons  tunneling 
between  a  normal  metal  and  a  superconductor  through  a 
thin  insulating  barrier.  If  one  assumes  ( cf .  Bardeen^) 
that  the  matrix  element,  M,  that  appears  in  the  transition 
probability  for  electrons  tunneling  through  the  barrier 
is  constant  for  small  energy  differences,  then  the  reported 
humps  in  the  conductance  curves  indicate  structure  in  the 
density  of  states  distribution.  Solutions  of  the  integral 
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equation  which  gives  the  energy  gap  parameter.,  A,  as  a 

function  of  the  excitation  energy  E  have  been  carried. 

4  R 

out  by  Swihart  and  Culler.,  Fried.,  Huff  and  Schrieffer  . 

Both  of  these  predict  structure  in  the  density  of  states., 

in  particular  two  humps  occuring  at  energies  near  k0^ 

and  2  k0j),  where  0D  is  the  Debye  temperature.  However,, 

2 

Rowell.,  Chynoweth,  and  Phillips  observed  far  more 
extensive  harmonic  structure  in  the  density  of  states  of 
lead  than  had  been  theoretically  predicted.  They  observed 
this  structure  at  energies  A  +  nErp,  where  E^  is  the  energy 
of  the  transverse  acoustic  phonon  in  lead,  and  n  is  an 
integer.  Since  such  structure  has  until  now  only  been 
observed  in  lead  it  is  of  interest  to  examine  other  super¬ 
conductors.  As  one  expects  such  structure  to  be  most 
readily  observed  in  superconductors  having  a  low  0^,  the 
most  suitable  choices  appear  to  be  mercury,  indium,  and 
tin.  The  purpose  of  this  letter  is  to  describe  the  con¬ 
ductance  vs.  voltage  curve  for  indium  obtained  from  our 
tunneling  experiments  with  In-Al^O^-Al  junctions. 

The  relative  conductance,  (dl/dV)  /(dl/dV)nn, 
was  derived  from  measurements  at  temperatures  where  the 
indium,  but  not  the  aluminum  was  superconducting,  along 
with  similar  measurements  at  4°K  where  both  metals  were 
normal.  Our  method  involved  a  bridge  technique^  which 
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Is  very  sensitive  to  the  quantity  1  -  (dl/dV)  /(dl/dV)nn, 
the  latter  being  determined  to  within  ±  0.0001.  The  level 
of  the  relative  conductance  curve  with  respect  to  unity 
could  only  be  determined  to  an  accuracy  of  -  0.0005. 
Measurements  of  the  relative  conductance  at  energies  below 
about  6  x  10  J  ev  were  made  at  several  temperatures,  and 
agreed  with  those  of  Giaever,  Hart,  and  Megerle^ .  Our 
new  results  for  higher  energies  are  given  in  Fig.  1.  The 
specimens  were  prepared  in  a  manner  similar  to  that  des- 

7 

cribed  by  Giaever  ,  they  ranged  in  resistance  from  50  to 
230  ohms.  Thermal  contact  was  provided  by  helium  gas  at 
temperatures  above  the  lambda  point,  while  a  film  of 
superfluid  helium  formed  the  thermal  link  at  lower  temp¬ 
eratures.  The  relative  conductance  seemed  to  remain 
constant  and  equal  to  unity  at  voltages  larger  than  20  mv, 
this  however  could  not  be  accurately  ascertained  since 
heating  effects  start  to  become  noticeable  at  these  voltages. 
The  prominent  structure,  in  particular  the  minimum  at  V  14.5 
mv,  remained  constant  even  at  the  lowest  temperature  avail¬ 
able  to  us  (T  =  0.5°K),  which  was  well  below  the  transition 

temperature  of  aluminum.  This  implies  that  the  structure 

2 

is  not  of  the  type  observed  by  Rowell  et  al.  ,  since  there 
are  no  humps  and  there  is  no  change  in  the  energy  at  which 
the  minimum  appears  as  the  temperature  is  varied  ( ie .  the 
change  in  the  size  of  the  energy  gap  has  no  effect  on  the 
minimum) . 
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It  is  evident  from  the  scale  on  the  ordinate 

of  Pig.  1  that  If  any  multiphonon  effects  of  the  type 

2 

observed  in  lead  are  present  in  indium  then  they  must  be 
smaller  by  more  than  an  order  of  magnitude.  This  is  not 
altogether  unexpected  since  the  electron-phonon  inter¬ 
action  in  indium  is  far  weaker  than  in  lead.  The  observed 
structure  occurs  at  energies  above  the  Debye  energy 

_  o 

(k©^  ^  9,5  x  10” J  eV)  and  therefore  an  explanation  in 
terms  of  the  phonon  spectrum  seems  unlikely.  We  have  also 
considered  the  possibility  that  the  effect  may  have  been 
due  to  the  volume  change  associated  with  the  superconduct¬ 
ing  transition  on  the  basis  that  the  boundary  of  the  Brillouin 
zone  shifts  and  modifies  the  energy  contours  in  its  vicinity. 
This  is  again  an  unlikely  explanation  since  in  indium  the 
superconducting  phase  has  the  larger  volume.  It  appears 
that  a  better  theoretical  understanding  of  the  factors 
affecting  the  density  of  electron  states  in  a  superconductor 
than  is  presently  available  is  required  in  order  to  explain 
the  observed  effect  in  indium. 

The  authors  wish  to  express  their  gratitude  to 
Dr.  S.B.  Woods  for  his  constant  help  and  encouragement 
during  the  course  of  these  experiments. 
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FIG.  1.  Relative  conductance,  (dl/dV)ns/(dl/dV)nn, 
of  an  In-Al^O^-Al  tunnel  junction,  plotted  against  energy 
at  2°K.  Aq  is  the  limiting  value  of  the  energy  gap 
parameter  obtained  from  measurement  at  T  =  0.5°K.  The 
dashed  line  shows  the  ratio  of  superconducting  to  normal 
density  of  states  calculated  from  the  BCS  theory  for 
T  =  0°K. 


